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PHILIPPE DI FRANCESCO AND RINAT KEDEM 

Abstract. We consider the cluster algebra associated to the Q-system for Ar as a tool for relating Q- 
system solutions to all possible sets of initial data. Considered as a discrete integrable dynamical system, 
we show that the conserved quantities are partition functions of hard particles on certain weighted graphs 
determined by the choice of initial data. This allows us to interpret the solutions of the system as partition 
functions of Viennot's heaps on these graphs, or as partition functions of weighted paths on a dual graphs. 
The generating functions take the form of finite continued fractions. In this setting, the cluster mutations 
correspond to local rearrangements of the fractions which leave their final value unchanged. Finally, 
the general solutions of the Q-system are interpreted as partition functions for strongly non-intersecting 
families of lattice paths on target lattices. This expresses all cluster variables as manifestly positive Laurent 
polynomials of any initial data, thus proving the cluster positivity conjecture for the Ar Q-system. We 
also give the relation to domino tilings of deformed Aztec diamonds with defects. 
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1. Introduction 

The Ar Q-system is a recursion relation satisfied by characters of special irreducible finite-dimensional 
modules of the Lie algebra Ar |T5]. It is a discrete integrable dynamical system. On the other hand, the 
relations of the Q-system are mutations in a cluster algebra defined in [14] . One of the goals of this paper is 
to prove the positivity property of the corresponding cluster variables, by using the integrability property. 

We do this by solving the discrete integrable system, which can be mapped to several different types of 
statistical models: path models, heaps on graphs, or domino tilings. The choice initial conditions for the 
recursion relations determines the specific model, one for each choice of initial data. The Boltzmann weights 
are Laurent monomials in the initial data. Our construction gives an explicit solution of the Q-system as 
a function of any possible set of initial conditions. 

It is a conjecture of Fomin and Zelevinsky [9] that the cluster variables at any seed x of a cluster 
algebra, expressed as a function of the cluster variables of any other seed y of the algebra, are Laurent 
polynomials with non-negative coefficients. In the language of dynamical systems, the choice of seed 
variables y corresponds to the choice of initial conditions. 

The methods used in this paper appear to be new in the context of the positivity conjecture. They have 
the advantage that they give explicit solutions, and extend immediately to T-systems [6] and to certain 
integrable non-commutative cluster algebras introduced by Kontsevich [16l [7] . 

1.1. The Q-system. First, let us recall some definitions. Let Ir = {l,...,r}, and consider the family of 
commutative variables {Qa,n ■ a € Ir,n G Z}, related by a recursion relation of the form: 

(1-1) Qa,n + lQa,7i-l — Qa,n ~ Qa+l,nQa~l,7n Qo,n = Qr+l,n = 1, (n G Z, a G Ir) ■ 

A solution of this system is specified by giving a set of boundary conditions. 

The original Ar Q-system [15] is the recursion relation p.l[) (where n > 1), together with the boundary 
conditions 

(1.2) Qa,o = l, Qa,i = chV(a;„), a € Ir- 

Here, ^^(a;^) is one of the r fundamental representations of the Lie algebra slr+i. The solutions, written 
as functions of the initial variables, are the characters of the Ar Kirillov-Reshetikhin modules [15] : 

Qa,n = chV{nUJa). 

The boundary conditions (|1.2p are singular, in that Qa,-i — 0. All solutions are therefore polynomials 
in {Qa.i ■ OL G Ir} (Lemma 4.2 of j4]). In this paper, we relax this boundary condition. Moreover, as we 
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are interested in the positivity property of cluster algebras, we renormalize the variables in (jl.ip Fl Let 




In the rest of the paper we work only with the renormalized variables in (|1.3p . and henceforth we refer 
to Equation II. 41 as the Ar Q-system. 

1.2. Q-systems as Cluster algebras. Instead of specifying the boundary conditions of the form (|1.2I) . we 
may specify much more general boundary conditions by picking a set 2r variables in a consistent manner, 
and setting them to be formal variables. Any solution is then a function of these formal variables. In order 
to explain what we mean by "a consistent manner" , we can use the formulation of these recursion relations 
as mutations in a cluster algebra. 

It was shown in [Mj that the system (II. 4p can be expressed in terms of a cluster algebra 0: 

Theorem 1.1. |14j The equations (jl.4p for each a and n are mutations in the cluster algebra without 
coefficients, of rank 2r, defined by the seed (xo,i?o), where xq is a cluster variable and Bq is an exchange 
matrix, with 



where C is the Cartan matrix of Ar . 

A cluster algebra of rank n without coefhcients is defined as follows ,9j. Consider an an n-regular tree 
T„, with nodes connected via labeled edges. Each node is attached to n edges with distinct labels 1, ...,n. 
To each node t, we associate a seed consisting of a cluster variable and an exchange matrix. The cluster 
variable Xt has n components, Xt = {xtfi, ...,Xt^n), and the exchange matrix i?( is an n x n skew-symmetric 
matrix. 

Cluster variables at connected nodes are related by mutations given by the exchange matrix. Suppose 
node t is connected to node t' by an edge labeled k, then the seeds at these nodes are related by a mutation 
/ifc. The effect of this mutation is as follows: 



The cluster algebra is the commutative algebra over Q of the cluster variables. 

In the case of the cluster algebra defined in Theorem ll.il there is a subgraph Sr of T2r, which includes 
the node with seed (|1.5p , and is the maximal subgraph with the property that the mutation of the cluster 
variables along each edge of the graph is one of the equations (|1.4p . The union of the cluster variables over 
all nodes of the graph is the set {Ra,n ■ a G Ir,n gZ}^ 

To describe the seeds in Sr, note that the mutation of the variable Ra^n-i into uses the variables 

{Ra,n-i, Ra-i,n, Ra,n, Ra+i,n}- If thcsc are entries of the cluster variable Xt at node i, then it is possible 
to apply the mutation t — > t' which sends Ra,n-i Ra.n+i- With our indexing convention, this is the 
mutation /Zq if n is odd, ^a+r if n is even. Such a mutation will be called a "forward" mutation, as it 
increases the value of the index n — 1 — n + 1. (Cluster mutations are involutions, hence we also have 
"backward" mutations which send n + 1 ^ ?i — 1 . ) 



^Alternatively, we can introduce coefficients in the second term (see Appendix A in [4]). The two approaches are equivalent 
in this case. 

■^Note that this graph is different from the graphs described in 1141 [4], which were minimal rather than maximal. 



(1.5) 
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Figure 1 . The subgraph of Q-system mutations for A3 . The blue dot stands for the node 0. 

We see that any cluster variable in 9r has the properties: (i) The variable consists of pairs the form 
iRa,m^T Ra,m^+i), ct & Ir and (ii) rUa-i — 1 < TOq, < nia-i + 1 for all a E Ir- Define a path on the 
two-dimensional lattice in the (n, Q:)-plane, by connecting the points M — {(toq,, a)}^^]^. Then condition 
(ii) implies that M is a Motzkin path, with steps of type (1, 1), (—1, 1) and (0, 1). Therefore we have 

Lemma 1.2. The nodes of Sr are in bijection with Motzkin paths on the square lattice with r — 1 steps of 
the form (1, 1), (—1, 1) or (0, 1), connecting vertices with integer coordinates {n, a) such that 1 < a < r . 

We define xm to be the seed with the variables {Ra,ma+i : i = 0, 1, a € Ir}- Let Mq denote the Motzkin 
path with rua = 0, a € Ir- Then Xg = xmq- Forward mutations fia and ^.a+r act on Motzkin paths by 
increasing or decreasing one of the indices rua- If the resulting path is also a Motzkin path, such mutations 
are guaranteed to be of the type (11.41) . 

The graph Sr is an infinite strip, and we display a section of it for r = 3 in Figure [TJ where we have 
identified the nodes sharing the same cluster variables. Most of the results of this paper can be reduced to 
working with the fundamental domain of this strip. 

1.3. Outline of the paper. The Laurent property [TU] guarantees that the cluster variables are Laurent 
polynomials when expressed as functions of the cluster variables of any seed in the cluster algebra. It is a 
general conjecture [HI that the coefficients of the monomials in these Laurent polynomials are non-negative 
integers. In this paper, we prove positivity for the cluster variables in S^- 

The Q-system is a way of relating the solutions of (|1.4p to all possible sets of initial data. Mutations of 
the cluster algebra allow us to move within the set of possible initial data, as long as this initial data is 
of the form xm- Our aim is to give an explicit combinatorial description of Ra,n, for each choice of initial 
data Xm within this set, as the partition function for a statistical model with positive Boltzmann weights. 

We proceed using the following steps: 
Step 1. We give a combinatorial description of {_Ri.„, n G Z} as functions of the initial data Xq. All other 
variables {Ra,ni a > 1} are discrete Wronskians of 1—a < j < a— 1}. The variables {i?i,n, n € Z} 

satisfy a linear recursion relation with constant coefficients, which are the conserved quantities of the Q- 
system. Here, the integrability of the system plays a crucial role. 

The conserved quantities are partition functions of hard particles on a weighted graph G^, with weights 
which are Laurent monomials in the cluster variables Xq. Therefore the generating function F^^\t) = 
X]n>o ^^R^,n is equal to the partition function of Viennot's heaps on Gr- Thus, there is a simple expression 
for F[''\t) as a finite continued fraction. 

Elementary rearrangements of the continued fraction allow us to express f[^^ (t) as the generating func- 
tion for heaps on different graphs with different weights. Our goal is then to prove that these rearrangements 
are mutations of the initial data xm . 

Step 2. We use a heap-path correspondence to re-express {i?i,n(xo)}riez+ as partition functions for 
weighted paths on a dual graph Gr- 

Step 3. For each Motzkin path M, we construct a graph Fm and edge weights 2/e(M) which are Laurent 
monomials in the cluster variables xm- The variables {i?i,n}riez+ are expressed as the generating functions 
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for weighted paths on Fm- This proves the positivity conjecture for and the claim that continued 
fraction rearrangements are cluster mutations of initial data. 

Step 4. We use the discrete Wronskian expression for Ra^n to interpret these variables as partition 
functions for families of strongly non- intersecting paths on Fm- This is done by a generalization of the 
Lindstrom-Gessel-Viennot [19[ 111] determinant formula for the counting of non-intersecting lattice paths. 
This implies the positivity property for Ra^n as functions of xm- 

The paper is organized as follows. In Section [2J we show that the Q-system amounts to a discrete 
Wronskian equation for the i?i^„'s, and that with a > 1 are expressed as a x a discrete Wronskians 
of We deduce that i?i,„ satisfies a linear recursion relation with constant coefficients, interpreted as 

the conserved quantities of the Q-system. 

In Section|31 the conserved quantities are interpreted as partition functions for hard particles on a certain 
target graph Gr, with vertex weights which depend on the initial data Xq of (|1.5p . We rephrase the linear 
recursion relation for Ri^„ in terms of generating functions which are simple rational functions. 

In Section |4] we re-interpret these rational functions as partition functions for Viennot's heaps [24] on 
the same target graph Gr- This gives a first proof of the positivity of -Ri.„ as a Laurent polynomial of 
the initial data xq (jl.5|) . We also give an explicit expression for the generating function of i?i^„'s, as a 
finite continued fraction. We show how a simple rearrangement lemma for fractions allows, by iterative 
use, to rewrite the partition function for heaps on Gr as a partition function for heaps on other graphs, 
with weights accordingly transformed. The proof that these rearrangements correspond to mutations of 
the cluster variables appears in Section |6l 

In Section [5] we reformulate the heap partition function in terms of weighted paths on "dual" target 
graphs, the language of which is more amenable to mutations. For each Motzkin path, we construct a 
corresponding graph with edge weights which are monomials in xm- Partition functions for paths on these 
graphs are the cluster variables For completeness, we also explain how to construct the dual graphs 

for the heap models from these graphs. 

In Section |6] we prove that Ri_n is the partition function for weighted paths on a target graph Fm with 
weights determined by the initial data xm ■ Using the transfer matrix formulation, we prove the statement 
that fraction rearrangements correspond to mutations of the initial data. This leads to the main positivity 
theorem for the cluster variables in terms of any initial data xm- 

The extension of this result to Ra,m a > 1 is presented in Section [71 where Ra,n is interpreted a 
la Lindstrom-Gessel-Viennot as the partition function for families of strongly non-intersecting weighted 
paths on the same target graph as for This completes the statistical- mechanical interpretation of 

all the solutions of the Q-system, and proves the positivity of all the cluster variables of the subgraph 
corresponding to the Q-system, when expressed in terms of cluster variables at any of its nodes. 

In Section [8] we discuss the limiting case of and present various exact and asymptotic path 

enumeration results, which correspond to picking initial data xm with entries all equal to 1. 

Finally, since the Q-system is a limit of the T-system, aka the octahedron equation, in Section [9| we 
give the relation of our results to the known results on the octahedron equation. We show how to relate 
weighted domino tilings of the Aztec diamond to our weighted non-intersecting families of lattice paths, 
and to interpret the result of cluster mutations in that language as weighted tilings of suitably deformed 
Aztec diamonds, by means of dominos and also pairs of square "defects" . 

The appendices include explicit examples of our constructions for the case A^. 

Acknowledgements: We thank M. Bergvelt, C. Krattenthaler, A. Postnikov, H. Thomas, and particu- 
larly S. Fomin for many interesting discussions. We thank the organizers of the semester "Combinatorial 
Representation Theory" and the Mathematical Science Research Institute, Berkeley, CA, USA for hos- 
pitality, as well as the organizers of the program "Combinatorics and Statistical Physics" and the Erwin 
Schrodinger International Institute for Mathematical Physics, Vienna, Austria. R.K. acknowledges the hos- 
pitality of the Institut des Hautes Etudes Scientifiques, Bures-sur-Yvette, France. P. D.F. acknowledges 
the support of the ENIGMA research training network MRTN-CT-2004-5652, the ANR program GIMP, 
and the ESF program MISGAM. R. K.'s research is supported by NSF grants DMS 0500759 and 0802511. 



6 



PHILIPPE DI FRANCESCO AND RINAT KEDEM 



2. Properties of the Q-system 

2.1. The fundamental domain of seeds. Our goal is give explicit expressions for Ra.n as a function 
of any initial seed data xm- First, we use the symmetries of the Q-system to enable us to restrict our 
attention to a finite fundamental domain of initial seeds, parametrized by Motzkin paths which have a 
minimum node at 0. 

There are three obvious symmetries of the system. A symmetry a : {Ra.n} — ^ {Ra,7i} of (|1.4|) is a 
map with the property that Ra^n — /(x) then cr(i?Q,ri) ~ /(""(x)). Equation (|1.4[) is invariant under 
<y{Ra,n) = Ra-n+1- Therefore, 

Lemma 2.1. "Time reversal": 

if Ra.n — f{Rl,0, R2,0, Rr,0] Rl,l, R2,l, Rr,l), 

(2.1) then i?Q,^_„+i — i?2,i, i?r,i; ^i,o, ^2,0, ^r,o)- 
for all n lE Z and a (z Ir- 

We also have the reflection symmetry of (|1.4I) . a{Ra.n) = Rr+i-a,n- Therefore, 
Lemma 2.2. 

if Ra.n — f{Rl,0, R2,Q, Rr,0': RlS,R2A, Rrs) 

(2.2) then i?.r+i_c(,n = f{Rr,0, Rr~l,OT--, Rlfi] Rr,l, Rr-1,1, Rl,l)- 

for all n Zi and a E Ir. 

Finally, we have the translational invariance, <j{Ra,n) — Ra.n+k- 
Lemma 2.3. 

if Ra,n f{Rl,0, R2,0, Rr,0, Rl,l, R2,l, Rr,l), 

(2.3) then i?Q,_„+fc — f{Ri^k, R2,k, Rr,k', Rl,k+1, R2,k+1, Rr,k+l)- 

for all n, fc e Z and a £ Ir- 

This lemma may also be viewed as a special case of the substitution property of deformed Q-systems 
defined in |3]. 

For any Motzkin path M — (rrii, m^), let fj^i^x) denote the function of x such that 
Let M + fc = (mi + k, rrir + k). Then Equation (|2.3I) can be written as 

(2.4) Ra,n = fL^i^Mo) = f!^::!,i^M^.+k). 

More generally, 

(2-5) Ra,n+k = /l^+fc(xM) = /a^+fc(xMo) = /i^"'' (xMo+fc) = f^H^M+k)- 

Therefore, 

Theorem 2.4. Let Ra.n — fa^\'^M.), where f2^ is a positive Laurent polynomial ofx.-M, ct G Ir, n eZ. 
Then as a function o/xm+A;; k Cz Ij, Ra,n is also a positive Laurent polynomial, Ra,n — /a^n^'^'CxM+fc)- 
We can thus restrict our attention to fundamental domain: 

Definition 2.5. The fundamental domain S'r for the Ar Q -system is indexed by the Motzkin paths with 
r — 1 steps of the form {{ma, a)}aeir' with Mina£i^{ma) — 0. 

There are exactly 3''""'^ Motzkin paths in the fundamental domain. As an example. Figure [5] shows the 
9 Motzkin paths of the fundamental domain for A3, and the mutations relating them, c.f. the graph S3 in 
Figure [T] 

2.2. Discrete Wronskians. For any matrix M, let Mf^''"'^^ be the matrix obtained by it by removing 
rows ii, Tfc and columns ji, ...Tjm- 
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Figure 2. The Motzkin paths for the 9 seeds in the fundamental domain for the As case, and 
the mutations relating them. The leftmost vertices of each Motzkin path lie on the vertical axis 
x = 0. 

2.2.1. Pliicker relations. Let T be an n x (n + /j)-matrix. The Plucker relations for the minors of T are 

k 

P=l 

In particular, when k = 2, 

^2 7^ |2iai,a2| |2^bl,fc2| _ |2^bl,a2| |2iai,b2| _j_ |2^b2,a2 | \rpbi,ai 



Let n = r + 1, a2 = ji, &2 = 32, and {T)i^ai = (^i.ii, {T)i^bi 



Then Equation (g^l) 

gives the 



Desnanot-Jacobi formula for the minors of the matrix M = T"^^^^ of size (r + 1) x (r + 1): 
(2.8) |M| |MJ,^| = |M^| iMg^l - |M^| |M^|. 

2.2.2. Wronskian formula for i?i,„. Using ()1.4I) . it is possible to eliminate the variables {Ra,n}a>i in favor 
of The remaining equations determine in terms of the initial data. As a consequence, 

satisfies a linear recursion relation with constant coefficients. This can then be extended trivially to all 

Define the a x a matrix Ma^n with {Ma.n)i.j — Ri.n+i+j-i-a- That is, 

Ri ,n—a-\-2 Rl,n \ 

Rl,n-a+2 Rl,n-a+3 ' ' ' Rl,n+1 



(2.9) 



\ Rl,n Rl,n+1 ■■■ Rl.n+a-lJ 

and define the discrete Wronskian determinant to be Wa.n = |-^^Q,n|- 
Lemma 2.6. We have Ra.n — Wa,n- 

Proof. Applying the Desnanot-Jacobi formula (|2.8p to the (a + 1) x {a + 1) matrix M with entries 

~ Rl,n+i+j-a-2 

with the choice of rows ii = \,i2 — a + \ , and columns ji — 1, j2 = a + 1, we have 
(2.10) Wa,„+iW^„,„_i = „ + W„+i,„W'„_i,„ 

for any sequence i?i,n, and any a > 1, with the convention that Wo,n — 1 for all rt. The sequence Wa,n 
is the unique solution to eq. (|2.10l) such that Wo,,i = 1 and Wi_n — Ri,n- Comparing this to the Q-system 
(|1.4[) . we deduce that Ra.,i = Wa n, a ~ I, and the Lemma follows. □ 
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The boundary condition i?r+i,n 
Corollary 2.7. 



1 yields the foUowing polynomial relation for Ri „: 



(2.11) 



W, 



r+l,n 



Rl,n-r+l 



Rl,n 

Rl,7l+1 



= 1. 



Rl,n Rl,n+1 ■ ■ • Rl,n+r 

2.2.3. Integrals of motion. The determinant Wr+i^n is a discrete version of the Wronskian determinant 



In the theory of linear differential equations, the Wronskian of r linearly 



independent solutions to an r-th order linear differential equation is a constant. This is proved by differen- 
tiating the Wronskian and noting that a linear combination of its columns vanishes, due to the differential 
equation. Conversely, if the Wronskian is a (non-zero) constant (so that its columns are linearly indepen- 
dent), there exists a vanishing linear combination between the column vectors of its derivative, namely 



/j'-'"'' = '^j/i"' ^' 1 where the /'s are a linearly independent set of solutions of these equations. 

Theorem 2.8. The variables {Ri,n\nez satisfy a linear recursion relation involving r + 2 terms: 



(2.12) 



^ (-l)™c.r+i_,„i?i^„+„ = 0, n e Z, 



m— 



with the coefficients Cq = c^+i = 1, and with ci, C2, some constant (independent of n) coefficients 
determined by the initial conditions. 

Proof. In analogy with the continuous situation, consider the discrete derivative Wr+i,n+i — W^r-t-i.n ~ 0. 
Since Wr+i.n and Wr+i,ri+i have r identical columns {{Wr+i,n+i)i.j = {Wr+i^n)i+i,j, i £ Ir), 

(2.13) Wr+l,n+l — Wr+l.n — dct (-Rl.ri+j+j-r-1 — {—^Y Sj_r+lRl.7i+i-r-l) — 0. 

l<i j<r+l 

As a consequence, there exists a non-trivial linear combination of the columns of this difference which 
vanishes. From the form of the entries of these columns (in which the indices are shifted by —1 relative to 
each other) the coefficients of this linear combination are independent of n. □ 

3. Conserved quantities and Hard Particles 

3.1. Conserved quantities of the Q-system. Since Wr+i.n — 1, it is a conserved quantity, i.e. it is 
independent of n. More generally, we claim that there are r + I linearly independent conserved quantities, 
and therefore the Q-system is a discrete integrable system in the Liouville sense. 



Ci-l — \{Mr+2,n)r+2 'I 



Theorem 3.1. The following polynomials 

i = 0,...,r 

where Cr+i = cq = 1, are independent of n, and co,...,Cr are the linearly independent conserved quantities 
of the Ar Q-system. 

Proof. This follows from the fact that Wr+2.n = 0, as a consequence of the boundary condition Wr+i,n = 1 
and the Q-system relation. The conserved quantities are the minors of the expansion of this determinant 
with respect to the last row, as in Equation ()2.12p . We get only r -\- \ linearly independent minors, since 



Co = Wr+l,k-l = 1 = Wr+l^k = Cr+l- 

Example 3.2. For r = 1, we have 

Rl.,k-2 
Rl.k-l 



□ 



Cl 



Rl,k 
Rl,k+1 



Ri,i 
Ri,o 



1 



Ri,o 



Using the Q-system for Ai to eliminate Ri^k+2 — {Rf 
we get the conservation law: 

Ri,k , 1 , Rl,k~l 



Ri,oRi,i Ri.i 
l)/i?i,fc_i and Ri^k-2 



(^i,fc+i 



1)/^L 



k-1, 



Cl 



Rl.k-1 Rl.k-lRl,k 



_ RiA ^ 1 ^ ^1,0 
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This is a two-term recursion relation in k, whereas the Q-system is a three-term recursion. The former is 
an explicit discrete "first integral" of the latter. 



Another way of understanding the conserved quantities of Theorem l3.1l is via the translational invariance 
of the cluster algebra, expressed in Lemma [2.31 We get the following immediate 



Corollary 3.3. The quantities Ci, expressed in terms of the seerf xq, are conserved, namely: 

(3.1) Ci(i?i,o, • ■ . , i?r,o; ^1,1, • • ■ , Rr,l) — Ci{Ri^k, ■ ■ ■ , Rr.k', Rl,k+1, ■ ■ • , Rr,k+l) 

for all k (z Z and for i = 0, 1, 2, r. 

3.2. Recursion relations for discrete Wronskians with defects. We now derive explicit relations 
between the Ci and the initial data. It is useful to work in the context of the Q-system, which is 

obtained from the Ar system by relaxing the boundary condition i?r+i.n = 1: 

(3.2) ra,n+ira,n~l ^ rl^^n + 1^a+l,nra-l,n, ro,„ = 1 , 71 G Z, a > 1 . 

Again, as in Lemma 12.61 we have 

ra,n ■= det iri,n+i+j-l-a)- 

1<2,J<Q 

Define the a x a Wronskians with a "defect" at position a — m {Q < m < a, n E Z): 

(^^•3) Ca,m,n = ^^'^^^^ ''l,n+i+s„_„(i)-Q-l j Sm{j) ~ -|- 1 if j > TO 

where Co,o,n = 1 for all n. Note that Ca,o.n — fa.n and Ca.a,n — 7'a,n+i- In addition, q = Cr+i^i,k if 
'fa,n = Ra.m that is, when we impose the condition rr+i,„ = 1. 

Lemma 3.4. The Ca^m.n satisfy the following recursion relation: 

(3.4) 7"q— l,n7'a — l.n+lCa.m.n — ^a.n^a—l,n-\-lCa—l,m.n ~^ ^a—l,n'^a,7i+l^a — l,7n — l.n "t" 7"q, 717*0,714-1^0 — 2, m—l,n- 

Proof. Applying p.Sp to the ax a matrix M with entries Mij = ri^n+i+j-a, hj G {1, • ■ • , ct}, with ii — 1, 
*2 — Oi, ji = a ~ 771 and j2 — a, 

(3.5) 7"Q,,77-f-lCa — 2,m — l,n '^a—l^n-\-1^0L — l.m^n — Ta — l.n^OL—l^m.n-\-l- 

Using this to simplify the sum of the first and last terms on the r.h.s. of (|3.4p . must prove that 

(3.6) T a — l.n+l^a^m.n — T a^n+l^a—l^m—l^n ~t~ 7'a,nCa — l,m,n+l ■ 

Multiplying (|3.6p by ra^n+i, and using (|3.2p . 

where we have used again (13. 5p . with a — )• a + 1, to simplify the second line. 

The last equation follows from (|2.7p . using the a x a + 2 matrix T with entries T^.i = 5i,a and Tij = 
Ti,n+i+j-a-i (2 < J < « + 2 and 1 < i < a), with oi = 1, 02 = 2, 61 = a + 2 — jn and 62 = a + 2. Then 
eq. (|2.7p becomes 

^0,71+2^0 — 1,771 — 1,77 — 777 ti^iTq- 1 7i_|_i 1" a .n+l^a — l ,m ,n-\-l 

which is the desired relation. □ 



Define 

1^0,77 



a = 1, 2, ... 77, G Z. 



Equation (|3.4p can be written as 

(3.7) Ca^jji^n — 7^a,TiCa — 1^777^77 ^ 7.^0,77+1^0-1,777-1,77 ~t~ 1'o,77 7^o,t7+1Co — 2,777 — 1,77 ■ 

Together with the initial conditions co,o,77 — 1 and ci,o,j7 = 7^1,77, ci,i^„ = 7^1,77+1, (|3.7p gives Ca^77i,„ as a 
polynomial in the variables {wfe.77, ffc.77+i}i</c<o, of total degree a. In particular, we have 

(3.8) Cq._o, 77 = 7;i^7i7;2, 77 ••• 7^0,77, Cq^q^ti = tii,77+i7;2, 77+1 ••• ?^'o, 77+1- 
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1 2 4 6 2r-2 2r 2r+1 





3 5 2r-1 

Figure 3. The graph Gr, with 2r + 1 vertices labeled i = 1,2, 2r + 1. 
Next we introduce the quantities which we call weights, for reasons which will become clear below: 

/o n\ ^a,n+l ^a— l,n^Q,n+l ^Q + l,7i+l ^a— l,n^a+l,n+l ^ 
(3.9) y2a-l,n = = ■ , V2a,n = = , a>l. 

We define 
(3.10) 



f l,nW2,n • • • Va,n 

Then p.7p becomes 

Theorem 3.5. T/ie quantities Ca.m,n of p.lOp satisfy 

(3.11) Ca^rn,n — ~t~ y2a — l,nC'Q-_i^77j_1^7^ ~l~ ^20 — 2,nC'cK_2,m— l,n; 

with the y 's as in eq. (|3.9p . 

Together with the initial condition Cafl,n — 1, this gives {Ca.m,n} as polynomials of homogeneous degree 
^ in yk,n^s, with 1 < A: < 2a — 1. 

Example 3.6. The first few C 's for a = 0, 1, 2, 3 read: 

Co,0,n — 1, Clfi.n = 1, Cl^l.n — yi.n, 

C2,0,n — 1, C2,l,n ~ 2/1, n + 2/2,71 + 2/3,71, C2,2,ri = 2/1, "2/3, nj 

C3,0,n = 1, C'a.l," = 2/l,n + 2/2,n + 2/3,n + 2/4,n + 2/5,n, 

^3,2 ,n — 2/l,^^2/3,?l2/5,n ■ 

We apply the above results to the conserved quantities of the Q-system of Theorem 13.11 We identify 
Ra,n = Ta,n by imposing the boundary condition rr+i.,i = 1 for all n. Then Wr+i.n = ^/fr.n and 

(3.12) Vi^nV2,n---Vr+l,n = 1- 

Therefore, 

(3.13) Ci = Cr+l,i,k ~ Cr+l,i,k ~ Cr+l,i,0- 

In particular, we recover cq = Cr+i = 1 from the explicit expressions for Cr+i.o,n and Cr+i^r+i.n of Equation 
(|3.8p . together with (|3.12p . We note that Cr+i,i.k — Cr+i,i,o are independent of k for all j — 0, 1, r + 1, 
in other words we have the conservation laws: Cj{{ya,k}) = ({2/0,0}) for all k G "Z. 

The identification p. 131) gives an expression for Ci in terms of the initial data Xq: By iterative use of 
the recursion relations of Theorem 13.51 for n = and 1 < a < r + 1, we get expressions for Cr+i,m,o as 
polynomials of homogeneous degree m in the weights 2/q,o, 1 < a < 2r + 1. These involve only the entries 

of Xq. 

3.3. Conserved quantities as hard particle partition functions. The recursion relations of the pre- 
vious section lead directly to an interpretation of the quantities Ci as partition functions of hard particles 
on a graph, with weights which depend only on xq. 

Let Gr be the graph of Figure |31 When r = 1, Gi is the chain with 3 vertices. To each vertex labeled i 
in the graph, we assign the weight yi. 

Let G be a graph with vertices labeled by the index set /. 

Definition 3.7. A hard particle configuration C on G is a subset of I such that i,jGG only if there is 
no edge connecting vertices i and j in G. 
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12 4 5 




Figure 4. The configurations of hard particles on G2, with, from left to right, 0, 1, 2 or 3 particles. 



Let HP(G') be the set of all hard particle configurations on G. 

If we assign a weight to each vertex i (z I, the weight of a configuration C is w{C) = Yiiec Vi- 
partition function of hard particles on G is 

(3.14) Z^^\{y,})^ ^(^)- 

CeHP(G) 

If we limit the summation to the set of configurations fixed cardinality j, we have the j-particle partition 
function Z^^\ 

In the particular case G = Gr of Figure [3l we have the partition function of j hard particles on Gr, 
denoted by := Zj ^ {yi, ...,t/2r+i). These satisfy recursion relations, coming from the structure of 

Gf . 

(G ) 

Theorem 3.8. The partition functions Zj ' satisfy the following recursion relations: 

(3.15) Zf = Zf-'^ + y2r+iZf_l-'^ + y2rZf:\-'^ 

Proof. Depending on the occupation numbers of the last two vertices, three situations may occur: 

(1) Vertices 2r and 2r + l are both vacant. This is a configuration of j hard particles on G^-i, obtained 
by erasing these two vertices and their adjacent edges. 

(2) Vertex 2r + 1 is occupied, and hence the vertex 2r is empty. Such configurations are those of j — 1 
hard particles on Gr-i- 

(3) The vertex 2r is occupied, and hence vertices 2r+l, 2r— 1 and 2r— 2 are empty. Such configurations 
are those of j — 1 hard particles on Gr-2 obtained by erasing the vertices 2r — 2, 2r — 1, 2r, 2r + 1 
and their incident edges. 

Each of these occupation states gives rise to one of the terms on the right hand side of the equation. □ 
Example 3.9. For r = 2, the hard particle model on G2 has the partition function 

(3.16) Z^^^\yi,y2,y3,y4,y5) = i + {yi + y2 + y3 + y4 + y5) + {ym + ym + y2y5 + y3y5 + yiy5) + {yiy3y5), 

where the various terms correspond to the configurations depicted in Fig^ 

Theorem 3.10. The j-th conserved quantity Cj of the Q-system for is equal to the partition function 

(G ) 

Zj ' (yi,n, 2/2r+i.n) forj hard particles on the graph Gr, with vertex weights yt = yi.„ defined in eq. (|3.9p . 
for i — 1,2, 2r + 1, and for any choice of n € 1j. 
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R2,Q 


l_ ^1,1 


RiflR2,i 


-^2,1 


1 ^1,0 


R2.I 


^1,0 


r ~" 

-R2, 0^1,1 


r ~ 

-Ri, 0^1,1 


-R2, 0-^2,1 


Rl,0 


|_ ^2,1 


R2,oRi,i 




1 ^2,0 


Ri,i 


R2,0 


r ~ 

-Rl, 0^2,1 


r " 

^2,0-^24 


^i,o-Ri,i 



Proof. The equations (I3.1ip and p.lSp are identical upon taking j/j = „. Moreover the initial conditions 

(G ) (G ) 

are also identical, as Zq = Cr+i.o,n = 1- We deduce that Zj ' — Cr+ij.n =^ Cj, independently of n. □ 

Corollary 3.11. The conserved quantities Cj can be expressed in terms xq as the partition functions for j 
hard particles on Gr , with vertex weights 

{ \ Ra—l oRa 1 Ra — 1 oRa+1 1 

(3.17) 2/2Q-1 = ^ 5 l<a<r+l, ?/2q = — 5-^-^ 1 < a < r. 

Example 3.12. In the case G2 of Examvle \3.9l we have (with j/ij/sj/s — I): Cq = C3 — 1, and 

(3.18) Cl : 

(3.19) C2 : 

The two integrals of motion of the A2 Q -system correspond to writing Ci and 02 with the substitutions 
Ra,i ~^ Ra,i+n-ij « = 0, 1, a = 1,2. These yield a system of recursion relations involving only indices n 
and n — \, as opposed to the original Q-system, which involves the indices n — l,n and n + I. 

3.4. Generating functions. 

3.4.1. A generating function for Ri^n- It is useful to introduce generating functions. Define 

(3.20) Ff'(t) = ^i?i,„r 

ji>0 

Theorem 3.13. We have the relation 

(3.21) Fi^\t) . ^^"^'"'l , ^±R,,,M-iy-W 

Proof. Consider the product of series ^ X^i^o ( X]jio(~l)"'cji-'J . Then all terms of order r + 1 or 

higher in t vanish, due to Theorem 12.81 We are left with the terms of order 0, 1, r, the term of order j 
being exactly {—lydj. □ 

Example 3.14. For r = 1, we have ci — + ^j^^ ^ + from Example \8.SX and do — Rifi, 
di — ciRifi — = "%77~; hence 



(3.22) )(0 = R., W.^^-.on ) 

yWi.o Jii^o-'i.i.i lii.i J 

3.4.2. Generating function and hard particles. 



Theorem 3.15. 

(3.23) F'f\t)^Ri 



{r)u^_ D Z'''^'-\0,-ty2,-ty3,:.,-ty2r+i) 



Zi°^'){-tyi,-ty2., -tys, -ty2r+i) 
with yi as in (|3.17p . 

Proof. In the expression p. 211) . the denominator is the partition function Z^^^\—tyi, —ty2r+i), accord- 
ing to Corollarv l3.11l The numerator of F^^\t) / Rifi is X]j=o(~*)"''^i/^i,Oi where 
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We proceed as for the Cj. First, we relax the condition that Rr+i^n — Ij hence work with the Va^n-, the 
solutions of p.2p . Define 

m 

(3.24) 7?„.„,„ = V(-l)"-'a.vn^^i^. 

Then djjrx^Q = Dr+ij.n = -Dr+i,j,0: independently of n when we impose the condition r^+Ln — 1, due to 
p.l3p . Substituting the recursion relations p. lip into this expression, we obtain an analogous recursion 
relation for Da,m,o- 

m 

Da,m,0 — * ' {y2a-lCa-l,i-l,0 + y2a-2Ca^2,i-l,0 + Ca~l,i,o) 

i=0 

= 2/2a-l-Da-l,,n-l,0 + 2/2Q-2£'a-2,m-l,0 + Da-l^m,Oj 

with Cq__i.o = and yk := yk,Q- The initial values of D for a = 1 are i^i^o.o = Ci,o.o = 1 and I?i,i.o = 
C'l.i.o — J/iC'i.o.o = 2/2 + 2/3- Both coincide with the values of Ci_o,o and Ci^i^o, respectively, when restricted 
to yi = 0. As the recursion relation for Da^m,o is identical to that for Ca.m,o\yi=o^ we deduce that 
fyi.o — Ca.m, o|yi=o for all CK, 771. This relation remains true after imposing the condition (13.121) . Therefore 

dj — Ri^Q Dr+l,j,0 = Rl.O C'r+l j",o(0, J/2, 2/2r+l) 

with the y's as in Corollarv l3.11l We deduce that the numerator of F^^^ (t) /i?i,o is equal to the denominator 
of (|3.2ip . restricted to the value yi — 0. □ 



3.4.3. Translational invariance. From the translational invariance property of Lemma 12.31 we may easily 
deduce an invariance property for the generating function F^^\t). Let us first write Fl^ \t) as an explicit 
expression F^^\t) = <l>((i?Q^o)Q=i; iRa,i)a=i',t) involving only the initial data xq. 

Theorem 3.16. The generating function $ satisfies the following translation property: 

k-l 

^{Rl,0, ■ ■ ■ , Rr,Q', Rl,l, ■ ■ ■ , Rr,l',t) — ^ Rl.nt"' + t''^{Rl,k, ■ ■ ■ , Rr,k', Rl,k+1, ■ ■ ■ , Rr,k+l',t) 

n=0 

for all k>0. 

Proof We write F^^\t) — J2t=o Ri,nt"' + t'' J2n>o Ri,n+kt^\ and apply eq.(|2.4p to aU i?i,„+fc in the second 
term. □ 

4. POSITIVITY: A HEAP INTERPRETATION 

We now have an expression for the generating function of i?i.„ {n > 0) in terms of the ratio of two par- 
tition functions of hard particles. Positivity of the terms Ri,m when expressed in terms of the fundamental 
cluster variables Xq, follows from a theorem relating this ratio to the partition function of heaps. This 
interpretation also allows us to find an explicit formula for the generating function of cluster variables. 

4.1. Heaps. Given a graph G, heaps on G are defined as follows (see [24j and the beautiful expository 
article [l8j). 

The graph G is represented in the xy plane in R"^, and we attach half- lines parallel to the positive z-axis, 
originating at each vertex of G. (See Figured]). The vertices of G are endowed with a partial ordering. 

On each half-line above vertex i, we stack an arbitrary number of discs of radius Ri and thickness a, 
with a hole at their center, so that they can freely slide along the half-lines (gravity points in the negative 
z-direction). The disc radii Ri are such that the distance between any pair of adjacent vertices i,j of G 
is < Ri + Rj, while the distance between any pair of non-adjacent vertices i,j is > Ri + Rj. Thus, the 
order in which the various discs are stacked matters only on neighboring half-lines (i.e. with connected 
projections on G), but not on distant ones. 

For a given stack of discs, its foreground is the set of discs that touch the a;y-plane. A stack is said to 
be admissible if it is empty or if its foreground is reduced to one disc, positioned at a vertex of smallest 
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order. (Such configurations are also called "pyramids" in the heap jargon [21]) We will call such admissible 
stacks of discs heaps on G. 

To each heap h on G, we associate a weight W{h) — Yideh^i'^)' where the product extends over all 
discs d of the heap, and where the weight W{d) = Zi \i d is stacked above vertex i. The partition function 
for heaps on G is 

(4.1) E ^(^) 

h heap on G 

Theorem 4.1. 

Z'^^>{-Zi,-Z2,-Z-i, . . .) 

Proof. This follows from the general theory of heaps [24] [18] . We write 

(4.3)$(G) (i^^i) ^(G) -Z3, . . .) ^ ^ M^(/i) w{c) = (0, -Z3, • • •) = E ^(^') 

(ft C) c' with 

vortex 1 empty 

where the sums extend over pairs {h, c) made of a heap h and a configuration c of hard particles on G, and 
configurations c' of hard particles on G such that the vertex 1 remains unoccupied. 

We define an involution Lp between pairs [h, c) which reverses the sign of W{h) 'w{c). Let c o ft, be the 
heap obtained by replacing the particles of c by discs, and by adding those discs on top of h. We define 
the background of any heap h to be the foreground of the heap obtained by flipping the conflguration 
upside-down. In particular, the background of c o ft contains c. Let d be the disc in the background of 
CO h that has the smallest vertex index i. Then if i is occupied in c, we form c' by removing the particle 
at i, and ft' by adding a disc on top of the i vertex. If i is unoccupied in c, then we form c' by adding a 
particle at i, and ft' by removing the top disc at i. Finally if ft is empty and the vertex 1 is unoccupied in 
c, then we leave the pair unchanged. These rules define an involution (^(ft, c) = (ft',c'). By construction, 
we have W{h)w{c) = —W{h')w{c'), if {h',c') ^ {h,c). Hence the distinct pairs {{h,c),{h' ,c')) image of 
one-another under ip contribute zero to the sum on the l.h.s. of eq. (14.31) . and we are left only with the 
contribution of the fixed points of the involution. The latter correspond to the situation where h — $ and 
c has no particle at the vertex 1, producing the r.h.s. of (|4.3p . □ 



4.2. Positivity from heaps. Applying Theorem l4.1l to the case of G = Gr, and comparing the expressions 
p.23p and (|4.2I) . we arrive at the main theorem of this section, which allows to interpret the i?i_„ as partition 
functions for heaps. 

Theorem 4.2. The solution to the Q-system for n > is, up to a multiplicative factor Ri^, the 

partition function for configurations of heaps of n discs on Gr with weights Zi = tyi, with yi as in (|3.17|) . 
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depth 

- 1 - 



1 -tyi 



1 -tvz 



1 -ty3-ty4 



1 -tys-tvg 



2r-ai- 



l - t y2r-3~ * y2r- 



1 - 1 yzr.i- 1 y2r 



1 - ' l2x. 



Figure 6. The graph Gr, arranged into a hierarchical structure. The depth of the vertices is 
indicated on the center, and the corresponding terms of the continued fraction for the generating 
function F[^\t) are displayed on the right. 

Proof. Using Theorem I3.15[ we may rewrite F[^\t) / Ri^o exactly in the form of the r.h.s. of (|4.2p . with 

G = Gr and the weights Zi = tyi. We deduce that F^^\t) / Ri^o is the generating function for heaps on Gr 
with weight tyi per disc above the vertex i. The coefficient of i" in the corresponding series corresponds 
to heap configurations with exactly n discs. □ 

Corollary 4.3. For all n G Z, is a positive Laurent polynomial of the initial seed xq. 

Proof. By Theorem l4.2l i?i,n/i?i,o for rt > is a sum over heap configurations, each contributing a weight 
made of a product of y.^'s. The resulting sum is therefore a manifestly positive polynomial of the j/i's, 
themselves products of ratios of initial data. For n < 0, the result follows by applying the Lemma \2A\ □ 



4.3. Continued fraction expressions for generating functions. The heap interpretation allows us to 
write an explicit expression for F^^\t) as a rational function of t. 

Theorem 4.4. 

^1,0 



(4.4) 



Ft\t) 



1-t- 



VI 



l-tV3-t- 



1-'H2t-+1 



where yi are defined in (|3.17[) . 



Proof. Define a partial ordering on the vertices of Gr by their geodesic distance from vertex 1. Any 
nonempty heap on Gr is constructed by repeating the following two steps: 

(1) Stack one disc above vertex 1. 

(2) Construct a heap on the graph G'r — Gr\ {1} (the graph Gr without the vertex 1 and its incident 
edges) . 

Let c{t) be the generating function for heaps on GJ,, then clearly F^^\t) / Ri_o = 1/(1 — tyic{t)). Similarly, 
c{t) is found via a similar construction of heaps on GJ,. Clearly, c{t) = 1/(1 — ty2d{t)), where d{t) is the 
generating function for heaps on the graph G" = GJ, \ {2}. 

To find d{t), we note that G" has two minimal vertices, but they are connected, so that a heap on G" 
may have as its foreground either node 3 or node 4 but not both. Thus, d{t) = 1/(1 — tj/s — ty4e{t)), where 
e{t) is the generating function for heaps on G" \ {3, 4}. This procedure is iterated (see Figure^, resulting 
in Equation (|4.4p . □ 



When expanded as a power series in t, (|4.4p has manifestly positive coefficients which are Laurent 
polynomials in the y^'s. 
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• a' 


a + b 


= a' 


, b' 


a b 


= a' b' 




a c 


= a' c' 



Figure 7. The local transformation of Lemma 14.61 

4.4. Continued fraction rearrangements. One can rewrite the continued fraction expression for F^^ ■* (t) 
in various ways using two simple Lemmas: 

Lemma 4.5. For all a, b there is an identity of power series oft: 

1 a 



(4.5) 



= 1 + t 



1-ta- tb' 



and 



Lemma 4.6. For all a, b, c, d such that c and a + b ^ 0, 

, . ha' 
(4.6) 

where 
(4.7) 



1 — t ^ 1 — f " 



a' = a + b, b' 



be 
a + b' 



a'b' 



a + b 



c^b' + c' 



(4.8) 



b' + d' ^ , ^ 

For example, applying Lemma l4.5l to the expression (j4.4p with a — y\ and b — yijiX — ty^ ■ ■ ■), we have 



F 



\t)^Rifl(^l+tj 



- tyi - 



ty2 



where yi are defined in p.l7p . 
Example 4.7. For r = 1, we have 



(4.9) 



F^'\t) = R,,o{l 



t -Ri,i 
'-Ri,o 



*y2i-+i 



^1,0 



1-t- 



'1,1 

^1,0 



«1.0«1,1 



Figure [7] is a graphical interpretation of Lemma 14.61 It shows that the generating function of heaps on 
G can be rewritten as a generating function for heaps on G" with new weights. Note that a, b, c, d and the 
associated nodes might stand for composite generating functions and the associated subgraphs. 

Starting from the initial seed xq and its associated graph Gr^ repeated application of Lemmas 14.51 and 
14.61 on the expression (j4.8p produces new expressions for F^^'^ (t) with manifestly positive series expansions 
in t. 

Writing the expression in (j4.8p as F^^\t) = -Fmo(^o; ^)i let Xj\/ be another seed in the fundamen- 
tal domain. Then we claim that there is a sequence A of applications Lemmas 14.51 and 14.61 such that 
A{F^^\t)) = -F'm(xm)- It turns out that we can generate in this way all the mutations of Xq within the 
fundamental domain 9^^- The proof of this statement appears in Section [5] Here, we illustrate this result 
with the example of A2 . The example of is given in Appendix \X\ 

Example 4.8. Consider the case of A2. We present all the rearrangements of F^ (t) which correspond 
to mutations of the weights within the fundamental domain. 
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Figure 8. The graph heap formulation of the apphcation of Lemma 14.61 on the generating 
function F^^\t) / Ri^ (a), and on the generating function "I>(t) such that F^^\t) = i?i,o + t-Ri,i$(t) 
(b). The circled vertices correspond to a, b, c in the Lemma [4.61 



(1) -Fmo(xo) with -X-o = (Rifi, R2,0i Ri,iR2,i): Equation \4--8\ is 
(4.10) (t) = i?i.o + ^ 



l-tys 



with 



/. ,,N -Ri.i R2.1 Ri,oR2.i Ri.o R2,o 

(4.11) yi = — , 2/2 = „ „ , ys = , 2/4 = „ p , y5 = 

^1,0 ^1,0-Ki4 -tt2, 0^1,1 ^2,0-^2,1 Jri2S 

(2) F^2(Mo)(M2(xo)), where ^J.2{^o) = (^1,0, ^2,2; -Ri, 1.^2,1). t^se Lemma IJ1>\ with a = y^, b = y^, 
c — ?/5 and d = 0. Then 

(4-12) fP(^)-i^i.o+ ^ "^^-^ ... = , 



1 - tei - - — ^il^ 1 - 



where the second expression follows from application of Lemma \4-5\ with a — Xi. Here 



2^3 = 



X4 



Ris 


X2 = 2/2 


_ -^2,1 


Ri,o 


RifiRi.i 


a + b 


= 2/3 + 2/4 


-Ri, 0(^2,1 + .^1,1 
R2,iRi,iR2,o 


be 


2/42/5 _ 




a' " 


2/3 + 2/4 


R\ lO' i?24-R2,2 


ac 


y^yb 




~^ ~ 


2/3 +2/4 


Ri^ia' R2.2 



R2,iRi,i 



(4.13) X5 = c 

by use of the Q-system. This is an expression for the generating function in terms of{Ri^o, R2,2] R1.1R2.1) 
/X2(xo). tt has a manifestly positive series expansion in t. 

Figure\^ (a) illustrates this transformation. The generating function is interpreted in terms of 
that for heaps on H2 ■' 

(2) Z''"^^0,-tX2,-tX3,-tX4,-tX5) 

^1 v^j = Ri,o- 



Z^^^'){ — tXi, —tX2, —tX3, —tX4, —tx^) 
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Figure 9. The graphs encoding the -Ri,n's for the case A2, and the corresponding mutations of 
cluster variables. 

(3) F^j(]v[o)(mi(™o)), where ^i(xo) = (i?i,2, -R2,o; -Ri,i^2,i)- Use Lemma [^7g| on ()4.10|) . with a = yi, 
b^y2,c^y3 + 2/4/(1 - ty,) = J^;;^^^ + i?i,i/(l - t^)) , and d ^ (see Figm (h)): 

(4.14) F[^\t) = i?i,o + t- 



1-t 



with 



Zl 



I , , , -^1.1 + ^2,1 i?i,2 

a = a + b = yi+y2 = 



1-tzr, a + b V l-ty^J R2,oRi,2Ri,i \ ' 1 - 

Z4 , ac f . \ -^1,1 / n2 , -^1,1 



^'•^'^ ^^ + T^ = ^ = ^ = i^'^ + ^ i?2,0i?2,li?l,2 ""^'^ ■ 1 - tP 

where we have used the Q-system. This yields positivity of all in terms o/(i?i.2, i?2,o; -Ri, 1^2,1) = 
/ii(xo), &2/ noting that Ri q = {R^ ^ + i?2,i)/-Ri,2- 

Figure \B[ (b) illustrates this transformation. (j4.14l) may be rewritten in terms of the generating 
function of heaps on K2 as: 

I?i2)^-^\ D , D + Z''^''^0,-tZ2,-tZ3,-tZ4,-tZ5,-tZQ) 

^1 V') — ^1 + ^1 I'' — fir w 

' Zy^2)(^-tZi^-tZ2,-tZ3,-tZi,-tZ5,-tZ6) 

Note that the graph K2 has one more node than G2 and H2, but the weights depend on the same 
number (A) of free parameters, since Z1Z3Z5 — 1, and Z2Z4 = z^zq. 

The rearranged expressions for F^ (t) considered above have been related to mutations of cluster variables 
in the fundamental domain J 2-: o,s well as to configurations of heaps on particular graphs, see Figure\^ 



5. Path generating functions 

In order to prove the claim of the last section, it is simpler to work with a path interpretation. There 
exist certain bijections between the partition function of heaps on a graph G, and the partition function of 
paths on an associated weighted rooted graph G. Such bijections are standard in the theory of heaps |24| . 

We will establish this bijection in the case of Gr- Then, we construct bijections between Motzkin paths 
representing cluster variables in the fundamental domain, and weighted graphs, on which the path partition 
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Figure 10. A heap on Gs with 8 discs, and the corresponding lattice path of length 16. The 
heap has small discs above vertices 3 and 5 of G3, and large discs above the other vertices. The 
large discs are in bijection with the descending steps, and the small discs are in bijection with the 
horizontal steps. This path has weight J/2j/3y42/5j/6j/7. On the left is the path graph G3, with 
labeled vertices (in black) and labeled edges (in gray). 

function gives the generating function Fl^^\t), in terms of the new cluster variables. For completeness, we 
give the bijection with the related graphs for heaps at the end of this section. 

5.1. Ftom heaps on Gr to paths on Gr- We start with the graph Gr defined in Figure |31 

Definition 5.1. The graph Gr is a vertical chain of r + 3 vertices labeled 0,l,2,...,r + 2, luith r + 2 vertical 
edges (i,i + 1) (0 < i < r + 1), together with r — 1 vertices 2\3' , ...,r' and r — 1 horizontal edges {i,i') 
(2 < i < r). It is rooted at its bottom vertex 0. (See the left of Fia \10\ for the r = 3 example.) 

The graph Gr is "dual" to Gr, in the sense that its edges are in bijection with the vertices of Gr- We 
may denote the edges by the same labels, i = 1,2, 2r + 1. 

Definition 5.2. The weights of Gr are defined as follows. There is a weight yi for each step along the 
edge i which goes towards the root vertex 0, and a weight 1 to all others. That is, the step i' i has weight 
y2i-i, the step i — > i — 1 has weight y2i z/ 1 < i < r + 1, step 1 — >■ has weight yi, and step r + 2 ^ r + 1 
has weight y2r+i ■ 

Paths along Gr, from the root to the root, also referred to as Gr-paths below, can be represented on 
a two-dimensional lattice as follows. Paths of length 2n start at the point (0, 0), end at the point (2n, 0), 
and cannot go below or above r + 2. They may contain steps of the type {i,j) (i + 1, j ± 1) (if j ± 1 
is in the range 0, r + 2), and steps {i,j) -> (i + 2,j) (2 < j < r). 

Theorem 5.3. The heaps on Gr with n discs are in bijection with the Gr-paths of length 2n, and their 
partition functions are equal, with weights as in Definition \5.2l 

Proof. Given a heap on the graph Gr, we associate a large disc to each vertex along the "backbone" 
{1, 2, 4, 2r — 2, 2r, 2r + 1}. We associate a small disc to the other vertices, {3, 5, 2r — 1}. The overlap 
between discs is given by the graph Gr , namely two discs overlap if and only if the corresponding vertices 
of Gr are connected via an edge. These discs are represented by bars in the plane in Figure [TOl 

We draw the bars on the faces of a tilted square lattice, in such a way that they are in bijection with the 
descending (for a large disc) or horizontal (for a small disc) steps of the two-dimensional representation of 
a Gr-path. In this picture, the descending steps of the path are the north-east edges of the square faces 
in which the large discs sit, while the horizontal steps are the horizontal diagonals of the square faces in 
which the small discs sit. 

There is a unique way of placing the discs with this constraint. We decompose the heap into successive 
shells. Each shell has faces containing the discs placed from the bottom to the top and from left to right, 
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5 



4 



2 



1 



Figure 11. The bijection between heaps on K2 and paths on the "dual" graph K2- The discs 
of the heap give rise to polygons, whose top right edges are exactly the descending steps of the 
path. We have indicated the weights Zi of the various discs, that are transferred to the descending 
steps of the path. We have indicated the vertex labels on K2, and the corresponding dual edge 
labels on K2- 



the spaces in-between being covered with only up steps. In this way, each large disc corresponds to a 
descending step and each small one to a double horizontal step of the path. The correspondence of weights 
is clear: descending steps receive the weights of the corresponding large discs, while the horizontal steps 
receive the weights of the corresponding small discs. 

Conversely, given a G^-path, we associate bijectively a large disc to each descending step and a small 
disc to each horizontal step. The resulting configuration is a heap over Gr, as the path always ends with 
a descending step corresponding to a disc over vertex 1 of G^. □ 

Corollary 5.4. For n > 0, the solution of (|1.4|) . expressed in terms of the variables xq, is i?i,o times 
the partition function of Gr -paths of length 2n, with weights as in Definition \5.2[ 

(r) 

The partition function for heaps on other the graphs corresponding to fraction rearrangements 
can also be written in terms of paths on graphs. For example, it is easy to see that the path graph 
corresponding to the heap graph Hr (the chain with 2r + 1 vertices) is a vertical chain with 2r + 2 vertices, 
with edge {i,i + 1) of Hr corresponding to vertex i + 1 in Hr- The edge weights for descending edges are 
the same as the corresponding vertex. Ascending edges have weight 1. 

Example 5.5. Consider the case of A2 (see Fig\^. We have a path formulation of the partition function 
in terms of paths on G2 as a function of xq . Paths on the graph H2 correspond to the generating function 
in terms of /i2(xo). The path formulation for the cluster variable /ii(xo), corresponding to heaps on K2, is 
the partition function for graphs on the graph K2 in Figure Note that the edge labeled 6 is an oriented 
edge (with weight zq). 

In general, it is quite complicated to work out the direct bijection between G and G. Instead, we 
introduce a bijection between Motzkin paths in the fundamental domain and path graphs, and a bijection 
between the same Motzkin paths and heap graphs. This establishes an identification between path and 
heap partition functions in the cases of interest. 

5.2. Motzkin paths and path graphs. 

5.2.1. The graphs Fm- For any Motzkin path M in the fundamental domain, we associate a graph Fm by 
(1) decomposing M into "strictly descending" pieces; (2) associating a subgraph to each descending piece 
and (3) gluing the subgraphs. 

(1) Decomposition of M: Each path M, consists of strictly descending pieces mi, m2, m^, where 
m^ = {(x, y), (x — 1, y + 1), ...}. (We consider a single vertex to be a descending piece). These are 
separated by p — 1 ascending steps of type (i) the step (1,1) in the plane, or (ii) the step (0,1) in 
the plane. 

(2) Graphs for m^: For each strictly descending piece m with k vertices, we define a F(fc) Fm as 
follows. When fc = 1, m is a single vertex, and F(l) is a chain of four vertices, as shown in the top 
left of Figure [H 
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Figure 12. The graphs T{k) corresponding to strictly descending Motzkin paths with k = 
1, 2, 3, 4 vertices. We have indicated the extra edge labels (in black) and the vertex labels (gray) 
along the vertical chain for the case when i = j — 0. 




Motzkin path Gluing 

Figure 13. The graphs rm;_i and F^. are glued together according to whether the pieces 
rtii-i and mi are separated by a step of type (i) or (ii). 

When m contains k > 2 vertices, the graph Fm — F(fc) consists of a "skeleton" Gk, plus extra 
oriented edges j^iifj — i>l,j<k + 2 and i > 0. There is a total of fc(fc — l)/2 extra oriented 
edges in F(fc), which we label by the vertices they connect. See Figure [T^ 

Each F(fc) has four distinguished vertices denoted by {t, t' , b' , b) = {k + 2,k + 1, 1, 0). 
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■1 



Figure 14. The graph Fmo corresponding to the vertical Motzkin path Mo with r vertices. 
The path is decomposed into k isolated vertices, each corresponding to a chain of type r(l), glued 
as indicated. The resulting graph is the tree Gr- 
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(■) 



Figure 15. The 8 trees T2r+2(ji, 12, is) for r — 4. We have indicated the values of the i's 
under each tree, and the labeling of the edges (gray) and vertices (black). 



(3) Gluing the subgraphs: Ascending steps of M are of type (i) (1,1) or (ii) (0,1). The graph Fm 
is obtained by gluing the graphs Fmi, i = 1,2,..., p as follows. Denote the edge {t',t) of F^; by 
{t^,ti), and so forth. We then make the following identifications of vertices and edges (See Figure 

ini): 

• Type (i): Identify ti_i with b'^ and t[_i with bi. The common edge is represented vertically. 

• Type (ii): Identify U-i with bi and t'^_i with b'^. The common edge is represented horizontally. 
We denote this gluing procedure by the symbol "|". Thus, 

(5.1) Fm = FmilFmal • • • iTnifc- 

Finally, all vertices are renumbered sequentially from bottom to top. 



Example 5.6. Consider the case of an ascending Motzkin path containing steps of type (i) and (ii) only. 
It decomposes into r isolated vertices, corresponding to Fm; = r(l)- Por the path Mq, with steps only of 
type (ii), the result is Fmq = Gr (see Fig\14^. 

More generally, Fm is obtained by gluing r chains F(l) vertically or horizontally. 

The path M is composed of p vertical chains of lengths li, ...,lp > 1, with steps of type (ii) only with 
Y^i=i — These are separated by p — 1 steps of type (i). Each vertical chain corresponds to a graph Gi-, 
which are glued according to rule (i). This results in a tree which has a vertical chain of length '2+'^^(li + l), 
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and consecutive sequences of — 1 horizontal edges, separated by pairs of vertices. The top two and bottom 
two vertices have no horizontal edge attached, (see Fia \15\ for r = A). 

The vertices ij to which the horizontal edges are attached I < ii < i2 < • • ■ < is < 2r— s where s = r — p, 
have the property that ii — \ is odd, and ij+i — ij is odd for all j < s. In fact, 

j j 

{ii,i2, . . . , = {0, 1, 2, . . . , 2r - s + 1} \ Uo<j<p{j + + ^ + ^^i' 

i=l i=l 

Definition 5.7. The trees obtained in the previous example are denoted by T'2r+2(*ii ■■■,is)- 

The 2r + 1 edges of such trees are ordered from bottom to top and labeled 1, 2, .., 2r + 1, including the 
horizontal edges, which have labels ii + 1, 12 + 2, ig + s. 

Noting that the sequence ja = {ia + a — l)/2 satisfies l<ji<j2<-''<js<r — 1 without further 
constraint, we see that there are exactly (J"^ ) such trees for r and s fixed, hence a total of 2''"^ when we 
sum over s = 0, 1, 2, r — 1. Note that when s — 0, T2r+2{) — Hr, "dual" of the chain of 2r + 1 vertices 
Hj. introduced in Section l5.ll 

This example is important, as it illustrates all 2^~^ possible skeleton trees any Fm in the fundamental 
domain can have: 

Definition 5.8. The .skeleton tree associated to M is the tree Fm', where M' is obtained from M by 
replacing all the .steps (—1, 1) by vertical steps (0, 1). In other words, Fm' is obtained from Fm by removing 
all its extra down-pointing edges. 

5.2.2. Weights on Fm- We will compute the partition function of paths on Fm, for which purpose, we 
assign a weight to each oriented edge of Fm. An unoriented edge, in this context, is considered to be a 
pair of edges oriented in opposite directions. The edge {v,v') is considered to be an ascending edge if the 
distance of v' from the vertex is greater than that of v. Otherwise, it is a descending edge. We assign 
weight 1 to all ascending edges, and weight 2/e(M) to each descending edge e. 

Edges are labeled as follows. Given Fm, consider the skeleton tree Fm' = T2r+2(*i, «2, is) of Definition 
15.71 Its edges have labels 1, 2, 2r + 1, which we retain for the graph Fm, and we assign to the descending 
skeleton edge i the weights yi (these are independent variables, not necessarily related to the weights 
encountered earlier). 

The extra descending edges of Fm are labeled by the pairs of vertices they connect, with 1 < 

j + l<i<2r — s (see Figure [T2|) . The weights ytj corresponding to down-pointing edges e = (z, j) with 
i — j > 1 can be expressed in terms of the skeleton weights. In view of the gluing procedure, we restrict 
our attention to strictly descending Motzkin paths m with k vertices (see FigfT2|). There are k{k — l)/2 
descending edges of type {i—>-j),2<j + l<i<k + l. Then 

(5.2) (?/2j,%+2,j, •■•,2/fe+ij) oc (y2j+i,y2j+2,yj+3j+i,---,2/fe+i,j+i) for j = l,2,...,fc- 1 

The proportionality is via overall non-vanishing scalar factors, so these relations allow to express all the 
weights ytj with i > j + 1 in terms of those of the skeleton tree, 2/1,2/2, 2/2fe+i- In fact, denote 

(5.3) yi+i,i-^y2t, and yi,t := yi' ,i ^ y2t-i- 
Then Equations ()5.2[) are equivalent to 

(5.4) y%,jym,t = yi,£yni,j for i > m and j > £ 
That is, 

ni—1 y-ri—l 
i^,yi+i.£ _ lh=jy2e 

[h=j+iye,i [le=j+iy2i-i 

These relations between edge weights will play a crucial role in Section [7] below, when we discuss the 
path interpretation of Ra.n- They will also become clear when we express the effect of cluster mutations 
on the weighted graphs. 
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t-1 t t+1 



Figure 16. The path-heap correspondence. To each descending step {t,i) — > (t + of any 
path on Fm, we associate the parallelogram indicated. The corresponding heap disc is a segment 
on the vertical diagonal of the parallelogram. To each double horizontal step we associate the 
triangle indicated. The corresponding heap disc is on the vertical median of the triangle. 

5.3. Prom Motzkin paths to heap graphs. For completeness, we give the bijection between the path 
graphs Fm and the heap graphs dual to them. This subsection is not necessary for further computations 
in this paper. We construct a heap graph Gm for each Motzkin path M via a generalized path-heap 
correspondence. 

Given a graph Fm, we represent paths on the graph from to on the two-dimensional lattice as before. 
We advance by one step in the "time" (x) direction for each step in the path, and record the height of the 
vertex visited in the vertical coordinate. 

Associate to each descending step {t,i) {t + in this picture the parallellogram with vertices 
{{t, i) , {t + 1, j){t, j — l){t — 1, i — 1)} , and represent a vertical segment {{t,x), x £ [j — 1 + e, i — e]} for e > 
sufficiently small (see FigflG)). For double horizontal steps (steps of the form {t—1, i) {t, i) — (t+l, i)), we 
draw the half-diamond with vertices {{t—1, i), {t+l, i), {t,i—l)}, and the segment {{t,x), x £ [i— l-fe, i— e]} 
for e > sufficiently small. These segments represent the discs of the heap. 

The heap graph Gm encodes the overlap between these segments, the vertices of Gm are in bijection 
with the descending steps on Fm, and the edges connect any pair of steps such that the associated discs 
cannot freely slide horizontally without touching each-other. Weights on the vertices are assigned according 
to edge weights. 

To construct Gm directly from M, we proceed as for Fm- We associate a graph Gm; to each strictly 
decending piece rrii of M. We then glue these pieces according to the type of separating step between 
them. 

Let m be a strictly descending Motzkin path with k vertices. Let Gm = G(fc). Its vertices are indexed 
by the descending steps on F(fc). The descending steps on F(fc) are indexed by {i,j) for i = j = 2,3, k, 
i = j -f 1 = 1, 2, k + 2 and fc-|-l>i>j-|-l>2, hence G{k) has a total of - 1 -I- fc -|- 2 + fc(fc - l)/2 = 
(fc + l)(fc + 2)/2 vertices. The descending steps fc + 2 ^ fc + 1 and 1-^-0 are singled out, and form the top 
and bottom vertices of G{k), denoted by t and b. 

The descending step {i,j) {i > j) overlaps with any descending step (1x1,1) {m > £) such that m or 
i e [j, i] (or both). The set of all these overlapping descending steps {i, j) — (m, £) forms the edges of G{k). 
We have represented the graphs G(fc) for fc = 1, 2, 3 in FigfTTl together with the path graphs F(fc) of FigfT2l 
for fc = 1,2,3. 

The graphs Gmi_i and G^m sue glued in two possible ways, according to whether nii 1 and mi are 
separated by a step (1,1) (type (i)) or (0, 1) (type (ii)) (see Fig fTS)) : 

• Type (i): The top vertex ti-i of Gmi_i is identified with the bottom vertex bi of Gmi- 

• Type (ii): The top vertex ti-i of Gmi_i is identified with the bottom vertex bi of Gmi, and all the 
vertices connected to are connected to all the vertices connected to bi, via additional edges. 

Example 5.9. Consider /ii(G3) of Fiaure \S6\ of Avvendix^M (top of second column from left). The corre- 
sponding Motzkin path is {(1, 1), (0, 2), (0, 3)}, and has two strictly descending pieces mi = {(1, 1), (0, 2)} 
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Figure 17. The heap graphs G{k) corresponding to the strictly descending Motzkin paths with 
fc = 1, 2, 3 vertices, together with the corresponding path target graphs r(fc). The vertices of G{k) 
are indexed by the descending steps on r(fc). 
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Figure 18. The heap graphs Gm-^^ and Cm; corresponding to the strictly descending Motzkin 
paths mi_i and rrii are glued in two possible ways, according to whether the Motzkin paths are 
separated by (i) a step (1, 1) or (ii) a step (0, 1). Both involve identifying the top vertex of Gin-_-^ 
with the bottom vertex of Gm- , but in the case (ii) all the vertices connected to the bottom vertex 
of Cm; must be connected to all the vertices connected to the top vertex of Gm-^-^. 
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and m2 — {(0,3)}. The graphs Gmi o.nd correspond to the vertices denoted {1,2,3,4,5,8} and 

{5,6,7} respectively in Fia \3b\ As they are glued according to the rule (ii), all the vertices connected to 5 
in Gmi, namely 4 and 8, are connected to all the vertices connected to 5 in Gmg ; namely 6. 

Remark 5.10. Let us define the overlap between descending edges ofT-^ as follows: the descending edge 
{hi) (i ^ j) overlaps the descending edge {m,t} (m > £) if and only if m or £ € [j,i] (or both). We 
may now bypass the above gluing procedure by directly associating to Fm the graph Gm is follows: (i) the 
vertices ofG-M oltc in bisection with the descending edges on Fm and (ii) the edges of Gm. connect all pairs 
of vertices such that the corresponding descending edges overlap. 

The bijection between paths on Fm and heaps on Gm follows as before from decomposing any heap into 
shells of successive foregrounds, and the fact that there is a unique way of arranging the discs (and their 
surrounding parallelograms) on the square lattice from bottom to top and left to right in each shell, and 
filling the spaces in-between with ascending steps of the path. 



Path partition functions give a combinatorial interpretation for {i?i.„}, expressed in terms of any cluster 
seed in the fundamental domain. For each Motzkin path M, F[^\t) = X]n>o^"-^i." ^ path partition 
function on a graph Fm with positive weights depending on xm. This results in a positivity theorem for 
all {i?i,„}„gz, as well as explicit expressions for the generating function. 

6.1. Path partition functions in terms of transfer matrices. 

6.1.1. Path transfer matrix. Let J'^''(a,6) be the set of paths on the graph Fm starting at vertex a and 
ending at vertex b with n descending steps. The partition function is 



This can be computed by use of the transfer matrix Tm, the weighted incidence matrix of Fm, with an 
additional factor t per descending step. Its rows and columns are indexed by the vertices of Fm (the 
ordering on vertices is such that i < i' < i with non- vanishing entries 



where / is the identity matrix. 

Let Fm{t) := 2'm(0,0). The entry ((/ - Tm)"^)o,o is computed by row-reduction, using upper unitri- 
angular row operations to obtain a lower triangular matrix, then taking the (0, 0) entry of the latter. To 
be precise, we iterate the following procedure: 

(1) Let t be the last row in the matrix A. Define the matrix A' with entries 



6. Path interpretation of and cluster mutations 




(6.2) 





Aj.tAtj 
At.t 



i < t. 



Then A'^ ^ = for all i < t. 
(2) Truncate A' by deleting its last row and column. 
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Repeat this until the resuh is a 1 x 1 matrix, which is l/(A~^)o,o- 

For our particular set of graphs, the reduction procedure has a graphical interpretation. It gives the 
partition function in terms of a pruned graph with the top part removed, and replaced by a single "loop" 
with a different weight, corresponding to the partition function on the pruned branch. 

Example 6.1. The transfer matrix of the graph K2 of Fiaure \12\ is 



(6.4) 





tyi 
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0/ 



with vertex order (0, 1, 2, 2', 3, 4). Note that ^3.1 
/I 



The reduction of A 








ty:i 

1 



IS 



1-t iiI-TKj-'))o, 

This is identical to the inverse of the factor of i-Ri.i in 
Zi ^Vi, i = 1,2,3,4,5 and zg = ya.i- 



the second term of eq. (|4.14p . upon identifying 



6.1.2. Block structure. To see how this works in general, we now describe the structure of Tm- It consists 
of blocks which are put together according to the gluing procedure. 

Recall the decomposition of M into strictly descending pieces mi. Let F^. be the graph with its 
bottom and top vertices and edges removed. Define T^ to be the transfer matrix of F^.. Let also 
F^^ and F^^^^ denote the bottom and top vertices, glued via vertical edges to F^^ and F^^. We write 
Tm = ^'ma I iTmi 1 1 ' ' ' I l-Tmp+i to denote the gluing procedure. (Gluing via the operation 1 1 consists of adding 
a vertical edge or a horizontal edge and vertex.) 

The matrix Tm is obtained by gluing the diagonal blocks Tr' according to whether the corresponding 
graph gluing is via a (i) vertical or (ii) horizontal edge (see Figure [T5)) : 

• Type (i): The blocks Tr' ^ and Tp' occupy successive diagonal blocks in the matrix Tm (vertex t[_T^ of 
F,-_;^ is followed by vertex 6- of Tr;). They are "glued" by the addition of two matrix elements, (TM)6'.,t'._^ = 
1 and (TM)t;_^,fc; = ?/ = ^VKA-i- 

• Type (ii): We place the two blocks Tr'_^ and Tp' in the matrix Tm in such a way that the last row 
and column of Tp' ^ (with index t^_-^^) coincide with the first row and column of Tp' (with index 6,^). We 
then insert a new row and column labeled bi following t^_l, with nonvanishing entries (Tm)^;,*' ^ = 1 and 

(TM)t;_j,6. = y = tyb,.t'._-^- 

6.2. Reduction and continued fractions. 

6.2.1. General reduction process. Now consider the reduction procedure performed on the matrix A = 
/-Tm. 

(1) Let {t',t) denote the last two row indices of A. The bottom right 2x2 submatrix of / — Tm has 

1 -v\ 

_^ ^ j with y = tyt^t' = ty-zr- 
replaces At'^t' = 1 1 — ?/• 



the form 



Reduction erases the row and column t of A, and 
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case (i) 





Tr 

'^i-l 
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case (ii) 




(a) 

motzkin 
path 




(b) 

graph 




(C) 

transfer matrix 



Figure 19. Construction of the transfer matrix Tm in the case of (i) vertical and (ii) horizontal 
gluing of the subgraphs r'^._^ and F^. . (a) Is the motzkin path, (b) is the graph gluing and (c) 
the gluing of the two diagonal blocks Tp/ ^ and Tp/ . We have shaded the graph in (b) and 
the corresponding matrix elements of the block Tp/ in (c). The added element y stands for the 
weight tyt^a in both cases (i) and (ii). 



(2) Inductive step: After reduction of all rows with index j > t'^, the resulting matrix Tp, ||...||r' 
differs from the original only in the element (i^,t'J, denoted by 7i, where 



(l-7.)-^ = ((/-rr,,,,||r;,,||...||r;,,)-^)^,^^^,^ 



We then reduce all indices j > t[_^ + 1. This results in the matrix A' . Its lower right 2x2 matrix, 
indexed by the vertices a, 6, has the form: 

Case (i): (^^^ ^ where y — tyb^ and (1 — c)^^ = ^{I — TY'.\\...\\r'^^J^^'^ ■ One more 

reduction step eliminates row and column b and gives the new entry Aa.a 1 — jzi^- 

Case (ii): ^ ^ . One more step in the reduction gives Aa.a f-s- 1 — y — c. 

In both cases, the net result is a modification of A in which the bottom right element is 1 —7^-1, 
where 

(1-7.-1)-^ = ((/-rp,.,_,,iirii|...||r;,,)-0,,_^,,_^ • 

The bottom right element of the reduced transfer matrix is 7i_i equal to (i) y/{l — c) or (ii) y + c. 

6.2.2. The case of a strictly descending Motzkin path. Consider the transfer matrix corresponding to r(fc)' 
(Figure [T^). For later use, we use a more general matrix T'{k), which has r(fc)'j,^-^ j.^^ = c instead of 0. 
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With vertex order (1, 2, 2', 3, 3', k, k' , k + 1), this matrix is 



(6.5) 



T'(fc) = 
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ty2 










ty 4,1 





tyk,i 





tyk+1,1 
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tyi 





ty4,2 





tyk,2 





tyk + 1,2 
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1 








ty5 


tye 





■ tyk,3 





tyk+1,3 











1 











■ 
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ty? ■ 


■ iyfe,4 





tyk + 1,4 
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ty2k-i 


ty2k 























■ 1 








V 




















• 1 





c 



Reduction of the last index in A = {I — T'{k)) results in 

yk+i,i 



in row i < fc — 1; 

V2k-2 ^ 2/2fc-2 + ™^ ^ ~ 1^ 

^kk t— >■ 1 — in row fc. 



whereas reduction of the next index changes AJ,^ = 1 — ty2k/(X — c) ^ \ — ty2k-i — ^fi^- 
Let 

(6.6) </5fe(y2,---,2/2fe;y3,i,y4,i,---,yfe+i,fc-i;c) = (^(/ - r'(fc))"^^^ ^ 

As a result of the two reduction steps, tpk is replaced by tpk-i, with suitable substitutions of variables. 



Lemma 6.2. The function tp^. = — T-^(^j^y) J is determined by the recursion relation 
</'fe({2/j}j=2; {yj..i}'2<i+i<j<k+i]c) 

= "Pk-liiy] + 0j,2k~2—. L=2 ; {yj.^ + ''j,k-. h<r+l<j<k; ty2k-l + t- ), 

i — c i — c i — c 



with 

(6.7) ¥'1(2/2; ;c) = 

Example 6.3. For k — 2, we have 



1 -ty2 
-1 1 - c 



1,1 



1-t^' 

1 — c 



/ \ I ys 1 tyi 

'/'2(y2, 1/3,2/4; 2/3,1; c) = </3i y2 + :r^; -^ys ' 



l-c' 



1 -c 



ty2- 



For fc = 3, we have: 



I ^ ( I ^4,2 , 2/4,1 , 

V3 (2/2, 2/3, 2/4, 2/5, 2/6; 2/3, 1,2/4,1, 2/4,2 ;c) = </32 2/2 , 2/3 , 2/4 + :;2/3,i + -; :;^2/5 



l-c' 



l-c' 



%6 

1 - c 



ty2+ 



'i',3.1+ l-c 



l-tya- 



ty4 + - 



Thus, is a finite, multiply branched continued fraction. Its power series expansion in t has coefficients 
which are polynomials in the weights with no n- negative integer coefficients. This expression is still valid if 
we relax the proportionality conditions (|5.2I) . but if we use these for j = 1, we see that 

where Wk does not depend on 2/2,2/3,2/4- 
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(a) (b) 

Figure 20. A path M with r = 9 (a), and the product of mutations (circled) yielding M when 
applied on Mq. In (b), the two generic moves (i) and (ii) corresponding to the actions of the 
mutation are indicated, together with their "boundary" version, (iii) and (iv) when a = 1, and 
(v) when Q = r. 

6.3. Mutations as fraction rearrangements. 

6.3.1. Mutations and Motzkin paths. One can associate a unique sequence of mutations to each Motzkin 
path M using the following procedure. The path M = {{ma, a)}a<£i^ is considered on the lattice, and each 
lattice point (x, a) with I < x < rua corresponds to a mutation /Ja i^J■a+r) if x is odd (even), respectively. 
The compound mutation fi such that /i(xo) — xm is the product of these mutations read from bottom to 
top, and from right to left. For example, in Figure [20] (a) , fj, — ^ii^i2Mi5MiM2M3M4M5M6/^7/^9- 

This restricted set of mutations acting on Mq yields any path M in the fundamental domain. We need to 
use only the two elementary moves shown in Figure [501 (b) (i)-(ii), and their "boundary" versions (iii)-(v). 
Without loss of generality we may therefore restrict ourselves to this subset of mutations. 

To see how this set of mutations acts on the partition function, we compare the graphs and weights 
(rM,y(M)) and (rM',y(M')), where M' = ^„(M) or fia+AM). 

6.3.2. Bulk mutations of type (a) and (h). Consider first the two bulk moves (i)-(ii) of Figure 1201 (b). Their 
effect on the corresponding target graphs is shown in Figure [21] (a) and (b). The path M is decomposed 
into a vertex labeled a, and bottom and top pieces mi and m2, corresponding to graphs r(l), Fi and r2 
respectively. 

In Figure [21] (a), the edge labels of the piece r(l) corresponding to the central isolated vertex a of M 
are identified, they are 2a — 1, 2q! and 2a + 1 We find that all the independent edge weights in Fm and Fm' 
are identical, except for the three weights y2a-ii y2a, ?/2q+i of Fm- These transform to 2/20-1; 2/201 2/20+1 
in Fm' . As we shall see below, the new weights correspond to the mutation of cluster variables. 

The matrix Am = 1 — ^Pm; softer the reduction of the block corresponding to F2, becomes: 
(6.9) 
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Motzkin 
paths 



Target 



graphs 2oc-i 





M' 



(a) 




Figure 21. The two bulk moves of Fig |20l (b) acting on the paths and graphs. The edges 
involved in the transformation are indicated. The mutation (a) involves the three edges of the 
center piece r(l), and is independent of the structure of mi and m2. Mutation (b) involves four 
edges, and depends on the length of the strictly descending subpath above vertex a+1, increasing 
its length by 1. On the graph it transforms the block r(fc — 1) into r(fc), creating k new descending 
edges. 

with (1 - c')"^ ^ [{I - TrJ^^ (&2 is the bottom vertex of r2). Three more iterations of reduction 

replace the bottom right element of / — Tr^ with 



(6.10) 




/ = 1- 


- ty2a-l 


ty2a 

^ 1- 






Similarly, reduction of / 


^ ^r^, yields (after 


reducin 


g the T2 part) 





















(6.11) / 
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■■ 


- 1 
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-1 


1-c' 


-ty'2a + l 






\ ••• 








-1 
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with c' as above. Three more reduction iterations result in the updated bottom right element of / — Tr^, 



(6.12) 



/' = 1 



ty'2c 
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The partition functions of paths on Fm and T^.' niust coincide, they are the same partition function 
expressed in terms of different variables. They coincide if and only if the weights 2/2q-1' 2/201 2/20+1 
such that / = /' (all the other weights are equal) . Application of the rearrangement Lemma I4.6[ with 
a = V2a-i, b = y2a, c = y2a+i and d = c', gives / = /' if and only if: 



3.13) 



2/2a-l = 2/2a-l +2/2a, 



y2ay2a+l 



y2o 



2/2o 



2/2Q+1 



2/2a-l2/2o!+l 



y2c 



y2o 



We interpret this transformation as the effect of the mutation fia or fia+r on the skeleton weights of 
Fm, resulting in a rearrangement of the continued fraction F]v[(t) into Fm'(/:). 

The mutation in Figure [21] (b) is slightly more subtle, because it depends on the length of the strictly 
descending subpath in m above a, whose length is increased by 1. All independent edge weights are 
unaffected by the mutation, except those of edges 2q! — 1, 2a + 2 (recall that the other edge weights are 
ratios of skeleton weights yi,?/2, •■•,2/2r+i)- 

The transfer matrix is 



(6.14) 



/ 



Tm. = 



V 

















ty2a-\ 








00 ■■• 01 


























hi,y2a+2('=-i)|r2 



where ry2Q+i,i/2c.+2 ~ 1) the circled graph in Fm, on the bottom of Fig |21l (b), for which we indicate 
the values of the weights of its first two lowest edges. After reduction, we are left with (/ — Tpi), with the 
bottom diagonal element replaced by 



(6.15) 



/ = 1 - ty2a-l 



ty2a 



1 - c(y2a+l, J/2a+2) 



where (1 - c(y2a+i, ^20+2)) ^ = [(^ ^ 

vertex of F(A: — 1). 
On the other hand, 



(fe_i)|r2) ) , where b' denotes the next-to-bottom 

/ b' .b' 



(6.16) 



7m' — 



V 


























■•■ 1 





ty'2^ * * ■ • • 


* 




1 













J/2Q + 1, 2/2a+2 ^ 


-i)|r2 



where the entries * stand for matrix elements of the form ty'j^^, where the vector (iy2Qj iWj.i)) proportional 
to the vector (^2/2q+ij ^2/2q+2j (^2/i,i)), which are the weights appearing in the row below it. This fact can 
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be used to eliminate the * entries using the row below. Then the reduction gives the matrix 



(6.17) 



/-Tt 



where (l - c{y'^^^-^,y'2a+2, 
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-1 



■0-1 
■ 

(/-Tr 






ty2a-l 







1 + 



V2c 
V'2c, + 1 
-1 



ciy'2a + l,y'2a + 2) ) 



(fe-l)|I 



^2a + l'''2a + 2 

reduction steps replace the bottom right element of (/ — TrJ by 



r2 ) ""^ I , with c as in (j6.15D . Two further 

/ h' ,h' 



(6.18) 



/' = 1- 



tyL-i 



1 



y2c 



2) 



The partition functions for graphs on Fm and Fm' must be equal, so we look for the transformation 

{y2a-l,y2a,y2a+l,y2a+2) ^ {yL-l ^ VL+l ^ y'2a+2) SUCh that / = /'. 

The function c' = 2(2/20+1,^2(1+2) identified by row-reducing (/ — Tj^(k)\r2)- Reducing the F2 part 



results in replacing the bottom right element of (/— rr(fe)) with 1 — d', where (1 — rf') 



This yields the relation 1 
Lemma 



V2c 



'Pk^iy2 



^2,^2 



2a,2/2Q + l,2/2a+2,2/2Q+3---,J/2a+2/c-2;{yj,i};d'), ^ith ifk aS in 

Therefore 



y'2c + l ~ 'Tk Vy2a7W2Q + l 

and we get c' = 1/20+1 + ''^fc 2/20+21 by comparison with the general expression 



(6.19) 



/ = 1 - ty2a-l - 



ty2a 



Applying the rearrangement Lemma 
that (|6.19p is satisfied if and only if 



1 - ty2a+l - ty2a+2Wk 

with a = y2a-i: 



/' = 1 



1 



Jl2a 



*!'2a + l+*y2c + 2™fc 



yL + l 1-*1'2c, + 1+*1'2q + 2'^'= 

y2a and c = y2a+i + 2/20+2 Wfe, we conclude 



(6.20) 



y2o-l = 2/2a-l + 2/2Q, 2/2a 



2/202/20+1 
?/2o-l + y2o 



2/20+1 



2/2q-12/2q + 1 
?/2o-l + ?/2o 



2/20+2 



2/2o+22/2o-l 
?/2o-l + y2o 



while all other weights y\ are equal to yi. This is interpreted as the effect of the mutation /io or /ir+o on 
the graph weights. 

6.3.3. Boundary mutations. We consider the mutations of FigHD] (b) (iii-v). 

• Case (v). We have a = r, which is just case (a) of Fig|?Tl with m.2 — and r2 reduced to a vertex. 
The transformation of weights is (I6.13p . 

• Cases (iii-iv). Here, a = 1. We take mi = and Fi a single vertex in both cases (a) and (b) of 
Figure 1211 However, the action of fii or fir+i introduces a new feature, which we call re-rooting. This is 
because the effect of the mutation is an application of the rearrangement Lemma 14.61 on the corresponding 
weighted path partition function FM.{t). This is possible only if the graph Fm is rooted at vertex b' = 1 

instead of vertex 6 = 0. We write Fm(<) = 1 + tyi(M)F!^{t), where FM(t) = ((I - Tm)^^) , while 

V / 0.0 



This re-rooting must take place whenever we act via the moves (iii-iv) as a direct consequence of the 
cases (a) and (b) of Figl^H Indeed, Fi is reduced to a vertex. The lower edge 2a — 1 = 1 is horizontal, and 
the vertex common with Fi is b' = 1, rather than 6 = 0, so that Fi = {b'}. The weight yi of the horizontal 
edge is associated with the step — )• 1 in the re-rooted formulation. 

In general, the weighted path partition function FM.{t), corresponding to the Motzkin path M, is related 
to the initial generating function F-M^^t) — R^QF^^\t), via the following sequence of re-rootings. We start 
from the canonical sequence of Motzkin paths leading from Mq to M = {mi, ...) via our restricted set 
of mutations. Within this sequence, the paths M'^**\ i — 1,2,..., mi are those which differ from their 
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predecessor only in via the action of the "boundary" mutations /ii or fir+i- Note that r7ii(Mf) — I. Only 
the boundary mutations [i\ and [ir+\ imply a re-rooting (otherwise -Fmc^-i) — -^m(»))- Thus, 

Lemma 6.4. The. 'partition junction F-^(i) is obtained from i^Mo(0 — ^iQ-^i^\i) sequence of 

re-rootings: 

FM,_At) = l + t2/i(M^-i)FM,W, 1,2,..., mi 
with i'M = -pMnij I where mi = mi(M) is t/ie first entry o/M = (m,Q, Q!)J^^i . 

This allows to interpret i?i.n+mi (xm), expressed as a function of the cluster variable xm, in terms of 
the partition function FM.{t). Let us denote by F{t) the coefficient of in the series F{t). 

Lemma 6.5. 

Ri.n+rn, = -Ri,oi^MoW = i?i.o yi (Mo)yi (Mi ) . . .yi (M„, J Fm (i) 

Proof. We use Lemma [631 The prefactor is obtained by collecting the successive multiplicative weights of 
each re-rooting. □ 

6.3.4. Main theorem. Let us summarize the results of this section in the following: 

Theorem 6.6. The mutation of cluster variables fia or fia+r ■ xm ^ xm' is equivalent to a rearrangement 
relating the continued fractions _Fm — >■ i^M' tho.t generate weighted paths on the rooted target graphs Fm and 
Fm'- The edge weights of the corresponding skeleton trees, y(M) — {j/i(M), ?/2r4-i(M)} andy' = y(M') 
are related through either of the two transformations ()6.13p or ()6.20p . while all other weights remain the 
same. 

6.4. Weights and the mutation matrix B. There is a simple expression for the edge weights of Fm 
in terms of the cluster variables xm and the mutation matrix _B(M) at the seed M. To specify all the 
edge weights Fm, one need only specify 2/i(M), y2(M), y2r+i(M), due to the relations ()5.2p for the other 
weights. 

Theorem 6.7. The weights j/i(M), y2(M), 2/2r+i(M) of the skeleton tree Fm are the following Laurent 
monomials of the cluster variable Xm- 

(6.21) y2a-i(M) = Z"'"" , (a = l,2,...,r-H), 

^5.22) y2a(M) = ( 1 - 0,n«.m^ + i + l(l 



X ( 1 - '^m„_i,m„ + l(l " ) ) , (a = l,2,...,r). 



where 



^ r,o\ \ _ Ra,n+1 Ra 

^.^oj .^a,n — , t^a,n 



Proof. Weights are determined by their initial values p.l7p at seed Mq and by the recursion relations 
()6.13p (in the case nia-i = TOq, < ma+i) and (|6.20p (in the case nia-i = rria = rria+i) for each mutation 
in our restricted set, all other weights being invariant. 

Using the Q-system (|1.4p . one checks that the relations (|6.13p and (16.201) are satisfied by the weights 
(|63T|) and ((6:221) . □ 

Note that eqs. ()6.2ip - ()6.22p involve only the data {Ra,m^ , Ra,m^^i)a=i, the cluster variables at the seed 
M. 

Example 6.8. Consider ascending Motzkin paths as in Examvle \5.6[ The weights from Theorem \6.7\ are 

(6.24) 2/2a-l(M) = , ?/2a(M) = , 

^rv — 1 _rrj ^ h^n m.^. 
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since only mutations of type (a) of the previous subsection are used. In the particular case of Is/l — Mg, 
Equations (j6.24p reduce to (j3.17p . 

In the case of the strictly ascending path Mmax = (0, 1, 1), with cluster variables {Ra.a-i, Ra,a)a<£ir, 

Ic ot;^ /l\/r \ Ra—l,a—2Ra,a , n Ra—l,o:—lRa+l,a 
(6.25) X2a-1 ■= y2a-l(M.rnax) = ^ , X2a := y2a[M.r,iax) = ^ ^ ■ 

The graph FM^a^ — Hr, the chain of 2r + 2 vertices. The partition function is 



(6.26) ^^M„_(t) = ^^f'(t)-E^"^i 



n>0 1 



Remark 6.9. Recalling the discrete Wronskian expression Ra^n = '^^^i<i,j<a{Ri,n+i+j-a-i), we may view 
the result (|6.25p - (|6.26p as a particular case of a theorem of Stieltjes [23] (see also which expresses the 
formal power series expansion F{X) = X]fc>o(~-'^)'^'^fc/'^'^^^ around oo for any sequence a^, k G Z>o, as 
the continued fraction 

F(X) = 



where - {^l? / {^l^l-i) , hk+i = {^DViKK+i), <^nd 

(6.27) det (a„+,+,_2) 

1<Z,J</C 

are Hankel determinants involving the sequence Uk. Indeed, writing X = —^/t and taking Ok = Ri,k, we 
easily identify the two continued fraction expressions F{X) andF[''\t) of eg. dOH)) upon taking Xk — p^p^^^ ' 
while = flQ = Ri^o yields the overall pref actor. Note that the continued fraction is actually finite, as 

The weights computed in Theorem 16.71 are simply related to the exchange matrix _B(M) of the cluster 
algebra at the seed M. 

Theorem 6.10. The exchange matrix _B(M) at the point M — {(ma, q:)}J^^x f^o,ds, for 1 < i,j < r: 

B(M),,,. = (-l)L^J ((-l)L^J - (-l)L'^J) 6\^-n,l 
B{M.\,,+r =(-l)L"^J + L^J+^a,j 
B(M),+,j = (-l)L¥J + L^Jc',j 

(6.28) B{m\+r,,+r = (-l)L^J ((-l)L^J - (-l)L^J) <5|,_,-|,i 
where C is the Cartan matrix of . 

Proof. By induction. The theorem is true when M = Mq, where -Bmo ~ \(j q 

Assuming it is true for some M with rUa-i — < rn^^i, let us prove it for M', with = mp + Sp^a- 
We distinguish according to the parity of TOq,: if is even, we use the mutation /i^, otherwise we use 

l^OL+r • 

Assume m„ even. Under the mutation the matrix elements _B(M)j ^ and _B(M)c, ^ are negated. 

This is compatible with Equation (lOS)) by noting that [^J = [^J and that = [^i^jtij + 1, 

which gives the expected extra minus sign. For the other entries of B{M), we use (|1.6p : 

S(M').,j = S(M),,, +sgn(B(M),,„)[B(M),,„S(M)„,,] + 

with i,j ^ a. Assume j < r. Then B{M)a,j ^ only if [^J = ^ ± 1, while |i - a| = \ j - a| = 1. Due 
to the Motzkin path condition, this is only possible if nij = — 1, but is impossible, as we must have 
ma+i,ma~i > ma. Therefore _Bjj(]V[) — Bij(M.'). Similarly, if i < r, we have the same conclusion. 
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If i,j > r, we write i ~ i' + r and j = j' + r, with i' ,j' < r. Then 

B{M),,+,,aB{M)^^^,+r = -C,.,„C„,,.(-1)L^J + L^J 

is positive only if i' = a and j' — ail, in which case (i) my = ma or (ii) j' ^ a + 1 and mji — 
rria+i ~ rria + 1, (or with i' and f interchanged). Then B{yV)i' j^r.aB{M.)a.j'+r = 2. In the case (i), 
we have Ba+r,j'+r{'^) = 0, hence Ba+r,j'+r{'^') =0 + 2 = 2, compatible with Equation (|6.28p . as 
l^m^j ^ Yi!^±L\^ +1. In the case (ii), we have Ba+r.a+i+r{^) = -2, hence Ba+r,a+i+r{^') = -2+2 = 0, 

also in agreement with Equation ()6.28p . as = ^ + 1 = [ "'"2"*"^ J • 

The case of ma odd is treated analogously. □ 

Remark 6.11. It is interesting to note that the B-matrices of Theorem ] 6. 1 0\ only have entries in {0, ±1, ±2}. 
This is true only for the cluster subgraph Sr, as the entries grow arbitrarily in the complete cluster graph 
T2r- We also note the remarkable property, that the sum of the four blocks r x r of the B-matrices always 
vanishes, namely that Bij + i?i+r,i + Bi_j+r + Bi^r.j+r = for i,j = 1, 2, r. 

The sequence (^{Ra,ma)a=iJ (^a,mo+i)a=i) ^-nd the sequence of cluster variables a;M = {Rcvcm Rodd) = 
{xi(M.))^Li (where the i?a,m are listed first for the even entries in m and then for the odd ones) are related 
via a permutation a-M- 

{a if ma is even 

a + r if ma is odd, and a < r 
a — r if ma is odd, and a > r 

We have the following expression: 

Theorem 6.12. The lueights computed in Theorem \6.7\ are related to the exchange matrix of the cluster 
algebra B{M): 

(6.30) ^Hf^^ = nx.(M)^(^)-."M(=), J^2^ = nx,(M)^(^)^.^M<"-), {a e /..) 

y2a{M) fj^ y2a+l{M) f-J^ 

6.5. Positivity of the cluster variables The variables i?i^„ can now be shown to be positive 

Laurent polynomials of any possible set of cluster variables xm- Note that the graphs Em are invariant 
under translations of M, but the same is not true the weights y(M). 

Theorem 6.13. For any Motzkin path M = {{ma,ot)\^a=i with r vertices and any n>Q, the solution 
^i,n+mi (xm) of the Ar Q-systcm, expressed as a function o/xm, is equal to Ri^mi times the generating 
function for weighted paths on Tm., starting and ending at the root, with n down steps, and with weights 
j/e(M) given by Theorem \6. 7[ namely: 



-Rl,ri+mi (xm) — Rl,mi ^m(0 



Moreover, Ri,n+mi is expressed as an explicit Laurent polynomial of the cluster variable xm, with non- 
negative integer coefficients, for all n G Z. 



Proof. By Theorem 12.41 we can restrict ourselves to the Motzkin paths of the fundamental domain. The 
first statement of the Theorem follows from Lemma 15751 with the prefactor i?i^o2/i(Mo)...j/i(Mp) = Ri.mn 
by use of (j6.2ip for a — 1. 

The second statement is clear for n > 0. It is a direct consequence of the reinterpretation of Ri^n+mi 
as the generating function for weighted paths with n down steps, for all n > 0, by also noting that the 
weights 2/e(M) are all positive Laurent monomials of the initial data, as a consequence of Theorem 16.71 

Eor < 0, we simply use an enhanced version of the reflection property of Lemma [2.11 Indeed, noting 
that the Q-system equations are invariant under the reflection n — )> —n, we deduce that the quantity 
Ri,n-i-mi, expressed as a function of the initial data xm, is the same as Ri^-n-mi expressed as a function of 
the reflected initial data Xm*, where the reflected Motzkin path M* = {(m^, a)} satisfies = — {ma + 1) 
for a = 1, 2, r. In other words, if Ri^n+mi — fi^M.), then i?i^_„_„i^ = /(xm*)- 
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Let us now use the translation invariance of the Q-system in the form of ea. (|2.5|) . We first translate 
the reflected Motzkin path M* back to the fundamental domain, namely so that its lowest entry be 

zero. This is done by considering to* = MinQ=i_2 rim-a), and the translate M' = M* — m*, with entries 

TO^ = TO^ — TO* for all a. Then for n < 0, the quantity Ri^^n-mi+m'^ is still given by the same expression 
as before in terms of the shifted Motzkin path M', namely i?i,_n-mi+m* = /(^m')- But now the first 
part of the Theorem applies, as to^ = toi(M') — — toi — 1 — to*. Indeed, as —n + 1 > 0, we deduce 
that Ri-n+i+m{ = /(xm') is a positive Laurent polynomial of the refiected-translated data xm'. As 
Ri.n+mi = /(xm) via the same function /, we conclude that Ri_n+mi is also a positive Laurent polynomial 
of the initial data xm, for all n < 0. This completes the proof of the Theorem. □ 

In view of the correspondence between path partition functions on Fm and heap partition functions on 
Gm, we have also 

Corollary 6.14. The weighted heap graph Gm associated to the Motzkin path M is constructed as above. 
The heaps on Gm with n discs are in bijection with the weighted paths on Fm with n descending steps, 
starting and ending at the root. For any Motzkin path M and any n > 0, the quantity Ri.n+mi is expressed 
in terms of the cluster variable xm as times the partition function for weighted heaps with n discs 

on Gm- 



We now provide a combinatorial interpretation of the determinant expressions for Ra^n (o^ > 1) as 
partition functions of families of strongly non-intersecting paths on Fm- For this, we need to generalize 
the usual notion of non-intersecting lattice paths. As a result we obtain a proof of the positivity of the 
determinant expression for Ra,n with a > 1. 

7.1. Ascending Motzkin paths: Paths on trees T2r+2{I)- First consider the case of the ascending 
Motzkin paths. The graphs Fm are skeleton trees, as in Section [5Jl denoted by T2r+2{I) where / = 
{ii, Z2, . . . , is] with 1 < ii < i2 < • • • < is < 2r — s and ii — 1, 12 — ii,. . - is — is-i odd. The variables 
are partition functions of paths on these trees, starting and ending at vertex 0. 

Paths on T2r+2 {I) are equivalent to paths on the two-dimensional lattice, with the y-coordinate restricted 
to0<?/<2r-|-l — s, with the following possible steps: (i) {x,y) — >■ {x + l,y+l); (ii) {x,y) — >■ {x+l,y—l) 
and (iii) {x, y) ^ {x + 2, y) (y e I). 

We say that two such paths intersect if they share a vertex. 

Theorem 7.1. Let n > a — I. Then Ra^n is (i?i_o)" times the partition function of families of a non- 
intersecting paths on T2r+2{I), starting at the points (0,0), (2,0), {2a — 2,0) and ending at the points 
(2n, 0), (2n -|- 2, 0), (2n -|- 2q; — 2, 0), with weights as in Examvle \6.8l 

Proof. We apply the Lindstrom-Gessel-Viennot (LGV) theorem [19l|Tl]. Consider the partition function 
■^Ai ' non- intersecting families of a weighted paths from the initial points Ai, A2, Aa to cndpoints 

Ei,E2, Ea. We assume that ii i < j and k < I, then a path Ai — >■ Ei must intersect any path Aj Ek- 
Let ZAi^Ej be the partition function for weighted paths from Ai to Ej, then 



Now, i?i,„ is Rifl times the partition function of paths on T2r+2{I) from (0,0) to (2n,0). Let Ai = 
{2i — 2, 0) and Ei — {2n + 2a — 2i, 0) {i — 1, 2, a). Since Za^^e, is the partition function for paths of 
2n -I- 2q; — 2i — 2j + 2 time-steps, it may be identified with Ri^n+a:-i-j+i/Ri,o- We conclude that 



7. Strongly non-intersecting path interpretation of R, 
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where we have used Lemma (|2.6p . 
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Figure 22. The eight pairs of non- intersecting Gr-paths of 8 and 4 time-steps for r > 3. The 
top box contains the non-intersecting Hr- (Dyck-) paths. There are six pairs of paths on G2. 

Example 7.2. For a = 2, n = 2 and r > 3, 

R2,3 = Rl,4:Rl.2 - (^1,3)^ = = xlx2X3{xiX3 + X1X4 + X2X4:) 

= yly2{yiyl + 2?/i2/3?/4 + yiyl + [vi + y2)yA{yb + ye)) 

with Xi as in (16.25^ and yi as in (j3.17l) . The first expression corresponds to the three pairs of paths in the 
top of Figure \2SX The second expression corresponds to the eight pairs of non-intersecting Gr-paths of 8 
and 4 time-steps depicted in Figure \2SX We also indicated in this figure the six pairs of non-intersecting 
paths on G2, for which no horizontal step at height 3 is allowed, eliminating the two configurations of the 
bottom row. 

7.2. Strongly non-intersecting paths on Fm- Theorem 17.11 can be generalized to paths on Fm even 
when Fm is not a tree. To do this, we must generalize the notion of non-intersecting paths, using a 
representation of Fjvi-paths on the two-dimensional lattice which preserves the property that the number 
of descending steps is half the horizontal length of the path. (Note that this representation is different 
from the one used in Figure 1161 ) 

7.2.1. Two-dimensional representation ofV-^-paths. Some of the descending steps on Fm may have height 
/i > 1, in case there is an edge j i with j ~ i > 1. When /i = or 1, each descending step is in bijection 
with an ascending one, but this is not true if /i > 1. We require that the total horizontal displacement in 
a path i— >i + l— j/— l—> 7— >i should be independent of the height h — j — i oi descending step 
used. Therefore, on the two-dimensional lattice, we draw a descent along an edge of length h as a, line from 
(x, y) to (x — h-\-2,y — h). Then a path which goes up h single steps, then down one step of height h has 
horizontal length 2, independently of h. 

With this convention, the horizontal distance between the start and end of a path from the to on 
Fm is twice the number of its descending steps (see Figure [23] for an illustration). 

Definition 7.3. The two-dimensional representation of a path with n descents on Fm is a path in (Z-|_)^, 
starting at (0,0) and ending at (2n, 0), with the possible steps: 

• {x, y) — > (x + 1, y + 1) whenever there is an edge ?/ — > y + 1 in Fm- 

• (x, y) ^ {x + 2, y) whenever there is an edge y ^ y' ^ y in Fm- 

• {x,y) ^ {x -\- 2 — h^y — h) , h > 1 whenever there is an edge y y — h in Fm- 

Note that these steps all preserve the parity of a; -|- y, which is even, so the path is on the even sublattice 
in Z^, with even x -\- y. 

7.2.2. Intersections of paths on the lattice. The determinant Ra,n = '^Q^i<i,j<a{Ri,n+a+i-i-j) may be 
interpreted using a generalization of the LGV formula. Intersections of the paths in Definition 17.31 are of 
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Figure 23. The two-dimensional path representation of a path on Fm with M = (2, 1,0) (see 
Figure [12] with r = 3). Descents of height 2 are vertical (horizontal displacement 2 — /i = 0), while 
descents of height 3 go back one step horizontally {2 — h = —1). With this choice, the horizontal 
distance between the starting and ending point is twice the number of descents (16 = 2x8 here). 



(h-2,h) 




h 



(0,0) 



Figure 24. The possible intersections of descending steps of two paths. The first descent is 
{h — 2, /i) — >■ (0, 0), and the second is of type (i + 2 — fc, i + fc) — >■ (i, i) with < k < h. There are 
h — k — 1 possible values of i. 



a more general type than the usual case. Paths may intersect along edges, and not just at a vertex. It is 
clear that edge intersections can occur only along descending edges. The possible crossing types can be 
catalogued as follows. 

Lemma 7.4. Consider a step of path Ti of type (h — 2,/i) — > (0,0). Then a descending step of a second 
path 5*2 niay cross this step only via the following h{h — l)/2 possible steps: 

(i + k-2,i + k) ^ {i,i), for < i < - fc - 1, 0<k<h-l. 

A generic case is represented in FiglMl 

7.2.3. A weight preserving involution on families of paths. Given the list of edge intersections in Lemma 
17.41 we define a weight-preserving involution on families of paths which we call flipping. Suppose two paths, 
Ti and ^2, intersect along an edge ei = {u + h — 2,v + h) — >• {u,v) of Ti and 62 = {u + i + k — 2,v + i + k) — >• 
{u + i,v + i) of J'2- Suppose = LiU CiU Ri (i ~ 1,2), where Li is the subpath of Ti before the vertex 
{u + h 2,v + h), and so forth. 
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Figure 25. A typical edge intersection of FM-paths (a) and the result of the flipping opera- 
tion on it (b). We have the weights of the steps. We have the identity of path weights, since 
yh+v,vyi+k+v,i+v ~ Vh+v ,i+vyi+k+v ,v The paths In (b) are said to be "too close" to each other. 

Definition 7.5. The flipping o/Ti and'J'2 w.r.t. the intersection along the edges ei and e2 the configuration 
of two new paths y'^, such that y[ = U U Ri, where 63 — {u + i + k — 2,v + i + k) — >■ {u,v) and 
e[ = {u + h — 2, V + h) ^ (u + i, V + i). 

Flipping is illustrated in Figure [55] The weight of two paths (Ti, J'2) is equal to the product of the 
weights of Ti and ^2 ■ 

Lemma 7.6. The weight of the two paths (Ti, J'2) is equal to that of the flipped paths (T'j^jCPj) of Definition 
Proof. 

w{ei)w{e2) = yh+v.vyi+k+v,i+v = yh+vA+vyi+k+v,v = w(e-)w(e2), 

By virtue of Equation (j5.4p . The rest of the weights of the path configurations are unchanged by the 
flipping operation, so the Lemma follows. □ 

For the graph Fm, there is a finite list of all possible results of the flipping procedure. 

Lemma 7.7. Given Fm; the results of flipping an intersection between a pair of paths is a pair of paths 
which include pairs of steps of the form ((u + h — 2,v + h) — > (u + i, v + i) , {u + j — 2, v + j) — > (m, v)) for 
any 1 < i < j < h, where there is a descending edge v + h v in Fm- 

Proof. By construction of Fm (see Section [5. 2. ip . if there exists a descending edge v + h v on Fm, then 
for all 1 < i < j < /i, the edges v + j ^ v + i exist on Fm, as well asv + h^v + i and t' + j — > v. The 
list of the definition therefore exhausts all possible cases of flippings of intersections along edges. □ 

Definition 7.8. Two paths obtained as the result of the flipping of an intersection are called "too close" 
to each other. 

For example, the paths in Figure [21] (b) are too close. 

Conversely, we may define the fiipping of a pair of paths which are too close to each other (with respect 
to the pair of edges that are too close) as the inverse of the transformation of Definition 17.51 The fiipping 
thus defined is an involution. 

Definition 7.9. Two paths are said to be strongly non-intersecting if (i) they do not intersect and (ii) they 
are not too close. 
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7.2.4. Generalization of LGV. We have the formula 

Ra,n+mi / Rl,n+mi+a+l—i—j \ 

{Rl,mi)°' ^<i.j<a \ Rl.mi J 

where n + mi > a — 1. Using Theorem 16. 131 and the path presentation above, Ri.n+mi+a+i-i-j / Ri,mi is 
the sum over weighted FM-paths from Ai = {2i — 2, 0) to Ej = (2n + 2a — 2j, 0). 

The proof of the LGV formula (j7.1l) uses an involution on configurations of paths, which leaves their 
weight invariant up to a sign. This cancels various pairs of intersecting paths in the determinant expan- 
sion. This involution is defined as follows: it leaves non-intersecting configurations invariant, otherwise, it 
interchanges the beginnings (until the first intersection) of the two leftmost paths that intersect first. The 
only terms not cancelled in the determinant correspond to the non-intersecting path configurations. 

The determinant (|7.3p is also as a sum over families of paths from {Aijf^j^ to the {Ej}'j^i, with their 
weight multiplied by the signature of the permutation a such that Ai is connected to E^(^iy 

We define an involution $ as follows. If a configuration has at least one intersection or edges which 
are too close to each other, then $ "flips" the beginnings of the two leftmost paths which intersect or 
are too-close, whichever comes first. That is, we either interchange the beginnings of paths if there is an 
intersection along a common vertex, or we apply the flipping procedure described above. Configurations 
which have no intersections nor edges which are too close to each other are invariant under $. 

The map <i> is clearly an involution. Applying Lemma 17.61 the total weight of a fiipped configuration 
is invariant. As for the ordinary case, $ pairs up terms in the expansion of the determinant which cancel 
each other, and we are left only with the ^-invariant configurations. Therefore the determinant ()7.3p is 
equal to the partition function of strongly non- intersecting paths of Definition 17.91 

Theorem 7.10. For any Motzkin path M, and for n + mi > a — 1, Ra.n+mi / {Ri,mi)'^ is the partition 
function of configurations of a strongly non-intersecting paths on Fm starting at the points Ai ~ {2i — 2,0), 
i = 1, 2, a and ending at the points Ej — {2n -t- 2q! — 2j, 0) , j — 1, 2, a. The weights are expressed in 
terms of the Motzkin path as in Theorem \6. 7\ and (jS.Sp . 

7.3. Positivity of Ra^n for all a, n. Theorem l7.10l implies that Ra,n+mi is a positive Laurent polynomial 
of any initial data, provided n -\- mi > a — 1, since all the weights in Theorem 16.71 and (15.51) are positive 
Laurent monomials of the initial data xm. 

We now consider the case when n + mi < a — 1. The determinant in eq. (j7.3l) involves some quantities 
Ri,m with indices to < 0, for which we have no interpretation as partition functions for paths. In order to 
generalize the result for arbitrary n € Z, we relate the expressions Ra,n+mi in the range n + mi < a — 1 
to some expressions Ra' ,n'+m' with n' + m'l > a' — 1. 

Consider the case when TOq, + 1 < n + mi < a — 1. Due to the structure of the Q-system, each Rp^n 
is inductively obtained from the values Rp^n-i, Ri3,n-2, -R^+i,n-i and i?^_i^„_i. Consequently, Ra^n+mi 
is a function only of the initial data that are contained in a "light-cone" of values of (/3,to), such that 
a — n — mi + k<l3<a + n-\- mi — k, for k < n + mi (sec Figl26p. 

In fact, Ra^n+mi dcpcnds only on the initial data xm' corresponding to a subset M' = {(ma, ol)}^^^^ 
of the Motzkin path M, namely such that (a, ma + 1) lies on the boundary of the light-cone, i.e. a = 
a — n — mi + k and TOq + 1 = fc. For the sake of this calculation we may freely modify the values of 
Ra-i,m and set them to 1, as they are not involved in the expression of Ra.n+mi- This has the effect of 
transforming the problem into one for Ar', with r' = r — a + 1. 

More precisely, Ra.n+mn as a function of a subset of the initial data xm, may be reinterpreted as the 
solution R'^, „,_|_„j/ of the Ar' Q-system, expressed in terms of the initial data xm', with = ma+p-i, 
/? = 1, 2, .., r — a + 1. In this new expression, we have n' + m'l = n + mi = a — a + ma -I- 1 = a' 4- TOq > a'. 

Note that the minimum to' of mp on the interval [a, r] may be strictly positive. In that case, we must use 
the translational invariance property of Lemma [2.31 (see also eq. (|2.4p ) to first translate both the Motzkin 
path M' and the index n' by —to'. We get n' — m' + m'l = a' + m'l — to' > a', and the Theorem 1 7 . 1 01 can be 
applied to conclude that R'^, n'-m'+m' ^ positive Laurent polynomial of the translated initial data xm" 
corresponding to to^" = m'^ — to' for all /3. By Lemma [2.31 we find that R'^, is a positive Laurent 

polynomial of the un-translatcd data xm' . 
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Figure 26. The generic situation rria + l<n<a— lin the (m, a) plane. We have represented 
the hght-cone of values determining Ra,n+mi (shaded area inside the wedge on the left of the 
point {a,n + mi)), as well as the initial data, in the form of two parallel Motzkin paths M = 
{{mp, j3)Yi3=\ and {{mp + l,l3)Yp^i- The light-cone cuts out a portion M' = {(rafi, l3)Yi3=a of M 
so that Ra,n+mi oxAj depcuds on the corresponding initial data. We may therefore truncate the 
picture by taking as new origin the (dashed) line (3 = a — 1, and interpret Ra,n+m.i as the solution 
Ra',n'+ma of ^hc Ar-a+1 Q-systcm, with initial data indexed by M', and with ct' = a — a + 1 and 
n' = n + mi — ma- 



We deduce that Ra,n is a positive Laurent polynomial of the initial data xm for all values of n > TOq, + 1. 

Finally, the positivity result is extended to n < (including negative values of n) by use of the 
corresponding invariance n — > — n of the Q-system (see Lemma |2.1|) . and the same trick as in the proof of 
Theorem 16. 131 involving the refiected-translated Motzkin path M'. We deduce the final: 

Theorem 7.11. The solution Ra.n of the Q-system for Ar is a Laurent polynomial with non-negative 
integer coefficients of any initial data xm indexed by any Motzkin path M, for all n e Z. 



8. ASYMPTOTICS 

In this section, we consider two limiting cases of the results: Solutions in the the limit r — > oo, corre- 
sponding to the Q-sytem of A00/2, and solutions when n 00. In the latter case, we are interested 
in the asymptotic behavior of the number of paths contributing to the partition function Ri_n as n — >■ 00. 

8.1. The limit Aao/2- In Equation (II. 4p . retaining only the boundary condition at a = 0, i?o,n ~ I {n £ Z) 
but dropping the condition at a — r + 1, gives us the solutions for the algebra Ar^^2- AH of the results of 
the previous sections generalize in a straightforward way. 

The initial seed is an infinite sequence, (i?Q,ma , ^Q.ma-i-i)aGZ+ indexed by semi-infinite Motzkin paths 
M. For example, for Mq — (0,0,...), and the corresponding heap graph is Goo/2: generalizing Figure [S] 
Then i?i^„ as a function of xmq is i?i,o times the generating function for weighted heaps on 0^0/2^ Fi{t) 
obtained as the limit of (|4.4p : 

^iW = ^ + ^ ^,. t 
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Using the rearrangement Lemmas 14.51 and 14.61 we can write the expressions for the generating function in 
terms of other initial seeds. For instance, the generating function corresponding to the "maximal" Motzkin 
path Mniax with rua = a — 1 is the continued fraction: 
1 



1-t 



yi 



i-t- 



The limit of the generating function corresponding to the strictly descending Motzkin path with r 
vertices, of the form 1/(1 — tyitpr), with as in Lemma l6.2l is a "continued fraction" with infinitely many 
branchings, a sort of self-similar object. 

8.2. Numbers of configurations. We consider the number of configurations contributing to i?i.„ in 
general. For this purpose, set i?a,m„ = Ra,ma,+i = 1 (o^ = 1; 2, ...,r). This implies that all the edge weights 
are j/e = 1- Therefore, Ra.n is a non- negative integer equal counts the numbers of configurations of the 
related statistical model. 

For example, when M = Mq, 

Lemma 8.1. For the initial data Ra.o = Ra,i — 1, a — l,2,...,r, the generating function Fl^\t) = 
I]n>o-^i>"^" reads: 

(8.3) Ft\t) - 1 + t^^^ ^ith P„(i) = t^C/„ (^-L - 

where U,n are the Chebyshev polynomials of the second kind, with 1/^(2 cos 0) = • "^^^ correspond- 

ing limit r — >■ oo for reads: 



(8.4) F^°"^ it) = l + tz{t) = 3 t VI 6t + t2 ^ ^ ^ ^ ^ + + 22t^ + + ZMt"" 



where z{t) — ^—^ — 6*+*^ . ^/jg generating function of the large Schroeder numbers. 

Proof. We use the expression l|4.4p with all weights equal to 1. Then there is a recursion relation, Pr+i{t) = 
(1 - t)Pr{t) - tPr-i{t), with Po{t) = 1 and Pi{t) = l-t. The limit z{t) := lim^^oo Pr{t) / Pr+i{t) therefore 
satisfies z = l — t+|. It also coincides with the continued fraction (18. ip with yi = 1 for alH. □ 

The rate of growth of considered as the number of paths of length 2n on Fmo: can also be 

analyzed. The radius of convergence of the series f[^^ (t) is given by the smallest root of the denominator 
of the fraction, which is Ur+i{^ — Vt). As the zeros of the Chebyshev polynomial Um are 2cos tt^^^^, 
fc = 1, 2, m, when n — >■ oo 

(8.5) i?i,„ ^ Cr X (^cos (-^) + ycos2(^)+l) 

for some constant Cr. 

The number of paths on T^.^^^ is also simple to compute, as it is the number of Dyck paths of length 
2n, which are limited by height 2r + 1: 



(8.6) F^^^Ht) 
In the limit r — )• oo. 



1 C/2 



(8.7) (t) = = l + t + 2t^ + 5t^ + Ut^ + 42t^ + 429t'^ 

is the generating function of the Catalan numbers c„ = (^j") . 
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As a result of the correspondence to domino tilings in the next section, this function counts the number 
of domino tilings of the (possibly truncated) halved Aztec diamond with d tiles missing. The function Ra.n 
counts the number of its indented versions. 

The large n behavior is i?i_„ ^ C'^(2cos{^:p^)) , for some constant C^. 

Now consider the expression as a function of the initial seed corresponding to M2, with = r ~ a 
(the maximal descending Motzkin path), where is as in Lemma[ 

1 



Flit) = + i'-'i?i,r-i<i>r(i;y), *,.(*; y) 



1 - iyi V'r ({yj J };c)' 

where (fr is defined in Lemma 16. 2 [ and the arguments are the weights j/i = ?/i(M2), i — 2,3, ...,2r, 
yij = yi,j{M2), 2 < j + 1 < i < r + 1, and c = ty2r+i(M2). 

Lemma 8.2. Setting Ra,r-a = Ra,r-a+i = 1 (a e Ir), 

(8.8) $,,(t)^l+<Ji(!) , Vrait) = i-irU^miVt). 

Vr+l[t) 

Proof. Setting the weights y's to 1 in 1/(1 — tipr), using the recursive definition of ipr of Lemma [6.21 with 
fm{{y]}-:{yi,]};ty2m+i)\y=i = i-e^{t) (™ - 0), = 1 and 6*0 = 0, we have Om = i!-t-e,^ii ■ Let 
1 — t — 6„i = ^y^^ , then we have a three-term recursion relation for Vm'. Kn+i = (2 — t)Vm — Vm-i, where 
Vo — l,Vi — 1 ~ t (or equivalently V-i = Vo = 1). The solution is Vm{t) = (-l)"'?/2m(v^), by iterating 
the Chebyshev recursion relation Um+i{x) = xUm{x) — Um-i{x), to get Un+i = {x^ — 2)C/„-i — Un-3, and 
identifying the initial conditions. □ 

The asymptotic behavior of with these boundary conditions is found from the smallest root of the 
denominator C/2r+2(v^)- Then 

(8.9) i?i,„ - c; X i ^ 



(2 Bin (2(2^ 



1-3). 

as n — )■ 00, for some constant C". 

In this case, the limit r — !• 00 of <J>r is ill-defined. The Motzkin path M2 has mi = r — 1 — cxd, hence 
does not have a good limit in our picture. More importantly, the graph has an infinite number of 
incoming edges at each vertex 1,2,3.... Hence, as we count paths according to their numbers n of down 
steps, we get an infinite number of paths from to as soon as n > 2. This can also be seen in the fact 
that the rate of growth in eq. (|8.9p diverges as 2r/7r. 

Example 8.3. In the cases r — 1,2,3, we have 

^ , , 1 l~3t + t^ UiiVi) 
*2(i = J = l + t - = l-t ^ ' 



" 1-t 



1 !t^— 1 - lot + 15*2 - 7*3 + t4 UsiVt) 



1-t 



More generally, for arbitrary choices of M and the associated initial conditions Ra,ma — Ra,ma+i — 1j 
we have F[^\t) = Lr{t) / Kr+i{t), where £„(<), Km(t) are polynomials of degree m with integer coefficients 
depending on M, and ^^(0) — 1 and Kra{t) ^ (^i)™ for large t. Indeed, from the results of Section [3J 
we know that Kr+i{t) is the generating function for hard particles on Gm, with weight ~t per particle. 
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So the empty configuration contributes Kr+i(0) = 1. Due to connectivity of Gm, the maximally occupied 
hard-particle configuration corresponds to all odd vertices 2« + l, i = 0,l,...,r being occupied by a particle 
(these are the duals of all the skeleton-tree edges of Fm, with odd labels). The corresponding weight is 
therefore {—tY'^^. 

For r large, it is always possible to reexpress Fi^\t) for any Motzkin path M obtained from Mq via 
finitely many mutations in the form of a rational fraction of the variables t and Pr-i(t)/ Pr{t). Indeed, the 
branchings of the continued fraction expression for ^ (t) all include some "tails" of the form 1/(1 — ij/a — 
t{yh/{l — tyc--.))) which, when we set all yt = 1, reduce to some ratio Pm{t) / Pm+i[t) . By use of the 3-term 
recursion relation for P„i, this can always be rewritten as a rational fraction of t and Pr-i/Pr- Then, in 
the case of ^00/2 > upon taking the limit r — 00, we see that, as Pr^i/Pr ^ z{t), F^^\t) takes the general 
form A{t) + B{t)^/l ~ 6t + t^, where A and B are rational fractions of t with integer coefficients. 

For instance let us take M — /iiMq. Then we have 



1 -t- 



1^3.1 

y2-^ TTfi — 



2 + t + 



where we have used the Q-system to rewrite i?i,o = ^ ' j^,^ ^ ' = 2, and the recursion relation for P to 
eliminate Pr-2- In the limit r 00, this yields for the Aao/2 Q-system with initial conditions Ri^i — 
Ri,2 — 1 and -Ra.o — Ra.i — 1 for all a > 2: 

„(oo).. „ , 1 - 5f + 2<2 + ^1 - 6t + i2 

Fl°^\t) ^ 2 + t + - — —^2 + t + t^ 



t - z{t) 1 - 7i + 5<2 _ t3 



9. The relation to domino tilings 
The solutions of the A^o T-system, also known as the octahedron equation [211 HOl EZ] , 

(9.1) Tij^k+iTij^k-i — Tij^i^kTij-i^k — Ti^ij^kTi-ij^k 

were given in j22j in terms of the partition function for domino tilings of the Aztec diamond and its 
generalizations. The Ar T-system (the same equation, with additional boundary conditions) is the fusion 
relation satisfied by the transfer matrix of the generalized, inhomogeneous Heisenberg spin chain jl5j . The 
Q-system (|l.ll) is a limit of the T-system, with the index j dropped, whereas the renormalized Q-system 
(|1.4p is obtained by formally dropping the index i. Therefore one should be able to recover the solutions 
of the Q-system from those of the T-system. 

Here, we give the explicit connection between the path formulation of Ra,n and domino tilings. In 
particular, we express Ra.ni^M.) for any Motzkin path M, as partition functions for tilings of certain 
domains of the square lattice by means of 2 x 1 and 1x2 dominos, and of rigid "defect" pairs of square 
tiles 1x1. 

9.1. Paths and matchings of the Aztec diamond. In this subsection, we consider the solutions Ra.n 
as functions of xq. For a restricted subset of the indices (a, n) there is an alternative combinatorial interpre- 
tation of Ra,n, related to the results of [22j on the solutions of the octahedron equation. This relation holds 
only for a restricted subset, because the boundary conditions of [22] are incompatible with our truncatiord 
Ro,n = Rr+i,n = 1- The conncction is therefore valid only sufficiently far away from these boundaries. 
Consider a system of equations with indices (a, n) €E Z x Z (the A^o Q-system): 

(9.2) Pa,n+lPa,n-l = Pa^n + Pa + l,nPa-l.n, «, n G Z 



'This truncation is also different from the truncation considered in 1121 for the so-called bounded octahedron recurrence. 
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To write Pa^n as a function of {pa.oj Pa.ijaez we use the recursion relation 

Pa.n+l 

Pa,n-l 

By induction, depends only on the subset of the initial data {p^OiP^.i} with a — 7i+l < /? < a + n~l. 
Hence, if we identify pa^i = Ra,i {i — 0, 1), then pa^n = Ra,n if 

(9.3) n < Min(a,r + 1 - a). 

We will show that Ra,m with indices in the set, is the generating function for positively weighted matchings 
of the Aztec diamond. 

Recall the interpretation of pa,7i in terms of partition functions of matchings of Aztec diamonds [22] . 
The Aztec diamond centered at (a, 0) of radius n is the set Aa,n = 7) S I |a — /?| + I7I < n} with 
boundary 2?q^„ = 7) g Z^, | |a — /3| + I7I = n} We denote by Aq_„ and -Bq,,„ the subsets of the square 
lattice with vertices in 1? + i) made of squares centered on the points of Aa,n and 23q^„ respectively. 

Consider the matchings, or compact dimer coverings, of A^ n- These are configurations of occupation 
of edges (including their vertices) of A^. ^ by dimers, where each vertex is covered by exactly one dimer. 
Each square (/3, 7) G A^^n has either 2,1, or edges covered, and we define e(/3, 7) = —1, 0, 1, respectively. 
Each square (/3, 7) € Ba,n has either or 1 edge covered, and we define e{fi, 7) = 0, 1 respectively. Let 

„ ,n\ fO if a + /3 + 7 + n = 0mod2 
^-"(^'^) = |l otherwise 

Then the generating function for matchings of Aa^n is 

(9-4) M^,n= n (p/3.e„.„(/3.7))' 

matchings of „ (/3,7)gyiQ,„uBa,„ 

The following is proved in |22) : 

Theorem 9.1. The solution pa^n of the system (j9.2p . expressed as a function 0/ {pct^Oi Pa.ijaez, is equal 
to the generating function for matchings of Aa,- 



Example 9.2. Figure \27\ shows, for n ~ 3, (a) the eight matchings of the Aztec diamond Aa^s, centered 
at (a,0) after rotation by 7r/4. The corresponding weights of e(j'. (|9.4p are (we denote by Oi the weight of 
the configuration labeled i in Fig\2^i 



_ Ra-2,lRa,lRa+2,l _ -^q-1,1-^q+1,1 _ 
"1 ~ R R ' "2 ~ p2 n ' "3 - ^ , 

-Rq-2,1-Rq + 1.1 -Rq-1,1-Rq+1,1 ^^a-l.l^a,l^Q+l.l 

04-7^ -p^ , 05-752"' "6-^ 5 5— , ay -as- -2 



^a-l,0^a+l,0 ^Q,0 -n-a-ljOJT-Q.O-flQ+i^o 



and we have Ra,3 = X]i=i'^2- Note that this expression is valid only for a = 2, 3,...,r — 1; r > 3, and 
provided we set i?o,n — Rr+i.n — 1 for all n. This can be compared to our previous result for Ra,3 as solution 
of the Ar Q-system, obtained by cutting out the initial data at P ~ a — i, hence working instead with Ar' , 
r' = r—a+3. This allows to reinterpret Ra,3 for Ar as i?3 3 for Ar-a+3, with initial data R'^, ^ = Ra'+a-3,i, 
i = 0, 1 and a' — 1, 2, r — a + 3. As such it is interpreted as {Ra-2,0)^ (= {R'l 0)^ ) times the generating 
function for triples of non-intersecting G^-paths (for any r > A) from the points (0,0), (2,0), (4,0) to the 
points (6,0), (8,0), (10,0), but with the weights y[ ^ bZ-To "■^'^ 2/2/3-1 r^XIZTorT+I'-Ti /"^ ^ 2 and 
2/2^ = fl|^^°~3'°flg^°~3'| fof /3 > 1. There are exactly 8 such triples, and their weights (multiplied by Ra-2 0) 
match one by one the weights Oi above. They are represented in Fig\21\ (b). 
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(b) 

Figure 27. (a) The 8 matchings of ^0,3 and (b) the non-intersecting Gr-paths from 
(0, 0), (2, 0), (4, 0) to (6, 0), (8, 0), (10, 0), where r > 4. The labeling corresponds to the weights at, 
i = 1,2,...,8. 




c d 

Figure 28. A bijection between lattice paths and domino tilings. The dual square lattice is 
bicolored, so there are 4 possible domino tiles, labeled a, b, c, d. Each corresponds to a step of 
path as indicated: a : (1, 1), b : (1, —1), c : no step, and d : (2, 0) — (1, 0) + (l, 0), double horizontal 
step. 

9.2. Lattice paths and domino tilings. There is a standard bijection between domino tilings and G,- 
lattice paths, (see e.g. [13]). Consider the dominos as tiling a chessboard. Then there are four types of 
2 X 1 or 1 X 2 dominos (see Figure [28]). These are in bijection with the four path steps: (1, 1), (1,-1), 
(2, 0) — (1, 0) + (1, 0) or no step at all. Thus, a domino tiling of a domain D can be rephrased in terms 
of a configuration of non-intersecting paths connecting points on the boundary dD. Conversely, consider 
Goo-lattice paths in the non-negative integer quadrant from (0, 0) to (2n, 0) of length 2n. (We call 
these Aztec paths in this section.) 

Lemma 9.3. Aztec paths are in bijection with domino tilings of the "halved" Aztec diamond HAn rep- 
resented in Figure f^P] of width 2n and height n + 1, with a floor at height h = —1/2, a ceiling at height 
h = n + 1/2, and a white face on the bottom left. 
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n+1 



2n-2 



Figure 29. The bijection between paths of length 2n steps from (0,0) to (2n, 0) and domino 
tihngs of a "halved" Aztec diamond HAn, of width 2n and height 




Figure 30. The local moves (i-ii) which give any Aztec path p from the maximal path po, and 
an example of such a transformation. First one goes from po to the Dyck path p closest to p by 
a sequence of "box removals" (i), before getting to p via "half-box removals" (ii). 

Proof. Any Aztec path may be obtained from the "maximal" Aztec path po up steps followed by n 
down steps) using the procedure shown in Figure 1301 

Given the path po, use the tiles a,b,d of Figure E51 to tile the region traversed by the path. One must 
then tile HA[^_i, the half Aztec diamond below this region, with width 2n — 2 and with the bottom row 
removed. The bottom left and right corner tiles are frozen, and are of type c, as are all the tiles touching 
the left and right boundaries of HA'^^_^. By induction, the entire domain HA[^_^ must be tiled with c 
tiles, a unique configuration. Therefore there is a unique tiling of HAn associated to the maximal path pq. 

For any other path p, we use the sequence of moves in Figure 1301 from po to p, and apply the corre- 
sponding transformations on the associated tiling by use of the bijection of Fig|28l This produces a unique 
tiling of HAn for each Aztec path p. 

Conversely, given any tiling of HAn, there is a unique path associated to it via the map in Figure B51 
A path can only enter the domain from its lower left end (1x2 domino of type a or 2 x 1 domino of type 
d), as all faces touching the left boundary are black, hence the tile d cannot be used. Similarly, the path 
can only exit via the lower right end (1x2 domino of type or 2 x 1 domino of type d). We deduce that 
there is a unique path associated to the tiling, that goes from (0,0) to {0,2n). □ 

We assign weights to the steps of Aztec paths according to the height h at which the step starts: 
• Weight 1 for steps of type (1, 1), irrespectively of h. 
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Figure 31. The bijection between families of a Aztec paths on Z+ from (0, 0), (2, 0), (2a — 
2, 0) to (2n, 0), (2n + 2, 0), 2n + 2a — 2, 0) and the domino tilings of the indented halved Aztec 
diamond IHAn,a of width 2n + 2a — 2 and height n + a. Here, a = 3. 



• Weight Z2h-2 for a step (1,-1) from height h to h ~ 1, h > 2; weight zi for a step (1,-1) from 
h = I to h ^ 0. 

• weight Z2h-i for a step (2,0) at height h, h > 2; weights wo,wi for steps (2,0) at heights h — 
and h — 1, respectively. 

Domino tilings receive weights according to the bijection of Figure [28l with domino c having weight 1. We 
will refer to those as weighted domino tilings in the following. 

9.3. Q-system solutions as functions of Xq and domain tilings. Recall the interpretation of i?i „ 
as the partition function of weighted paths on Gr- These are weighted Aztec paths with the following 
identification of weights: Z2r+i = wq — wi — Zj = (j > 2r+3), Zi = yi {i = 1, 2, 2r) and 22r+2 = 2/2r+i, 
with j/j as in (I3.17|) . The condition Zj = for j > 2r + 4 corresponds to truncating the tiled domain to the 
inside of a strip of height r + 3, with a ceiling at height h = r + 5/2, with floor at height h — —1/2. The 
conditions = wi — Z2r+i = Z2r+3 — forbid the use of the tile d in the two bottom (—1/2 < h < 3/2) 
and top (r + 1/2 < /i < r + 5/2) rows. 

Lemma 9.4. The families of a non-intersecting Aztec paths on from (0, 0), (2, 0), (2a — 2,0) to 
{2n, 0), {2n + 2, 0), 2n + 2a- 2, 0), with steps (1, 1), (1, -1) and double steps (2, 0) = (1, 0) + (1, 0), are 
in bijection with domino tilings of the "indented halved" Aztec diamond IHAn^a represented in Fia \31\ 

Proof. Along the same lines as for the a = \ case of Lemma 19.31 We first associate a tiling to any 
configuration of a non-intersecting Aztec paths. We simply use the fact that any configuration of a non- 
intersecting Aztec paths may be obtained from the configuration where all paths are maximal via successive 
applications of the local moves of FiglSOl We apply the construction of the proof of Lemma 19.31 to the 
bottom-most path, then to the next bottom-most, etc. 

Conversely, to construct the path configurations from the tilings, we note that the two successions of 
a — 1 indentations of the bottom boundary (left and right) impose the presence of a — 1 tiles a and one a 
or d on the left and a — 1 tiles b and one 6 or d on the right, and that paths can enter or exit nowhere else 
on the boundary, due to coloring constraints. This corresponds to exactly a paths touching the boundary 
of the domain, entirely determined by the tiling. □ 

Given that Ra^n is the partition function of a non-intersecting Aztec paths, with the identification of 
weights as above, this gives the relation to domino tilings. 

For sufficiently large r, the restriction of the weights effectively reduces the domain to be tiled to the 
shaded region in Figure 1311 as the tiling of the rest of the domain is entirely fixed (only dominos a on the 
left, and 6 on the right). The case a = 1 is also depicted with the same convention in Figl^ 

In the particular case when a = n, and r > 2n — 1, the shaded domain to be tiled is exactly the 
Aztec diamond itself (with width and height 2n ~ 2), with no constraints. It is interesting to compare 
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the partition function D„ {y) for the weighted tilings of this domain to the partition function for weighted 
Aztec diamond matchings (j9.4p . On one hand, we know that it is equal to 

(9.5) Dn{yi,y2,-.-,y2r+i) = -——^^^2 

where we have used the interpretation of Rn,n/Rio as generating function for G^-paths, applied the path- 
tiling bijection of Section [9?2l and removed the contributions of the tiles outside of the shaded domain. 

As discussed in Sections 17.31 and [9.H the solution pa,n of the Aoo Q-system for pairs a, n obeying ()9.3|) is 
identified with R'^ „, solution of the Ar Q-system for any r > 2n — 1, and with initial data R^ q = Ra-n+i3,o 
and R'j^ ^ = Ra-n+i3,i for /3 = 1, 2, r, namely with the weights (see also Example I9.2p : 

, Ri3+a-n,lRl3+a-n-l,0 r o i o i 1 

2/2/3-1 = ^ ^ for /3 1,2,...,?' + 1 

-ft^+a-n,0-ft/3+a-n-l,l 
, Rl3+a-n-l.oRp+a-n+l.l r a i o 

2/2/3 = — ^ ^ for/3 = l,2,...,r 

So we get two different ways of computing the same partition function p^.n for Aztec diamond matchings, 
one via the expression (|9.4p . the other as i?^ „. 

As is apparent from the Example 19.21 the weights leading to the expression 19 .41 and those of the families 
of non-intersecting paths are in bijection, but are incompatible with our weighted bijection. Indeed, if we 
rotate all the matchings of Figl27] (a) clockwise by a quarter-turn, we find that the corresponding dual 
domino tilings match the configurations of non-intersecting paths of Fig l27l (b) via our bijection only in 
the cases 2,3,4,7,8, while some permutation of 1,5,6 must be applied. This suggest perhaps that other 
natural bijections should exist. 

9.4. Solutions of the Q-system as functions of xm and tilings of domains with defects. We now 

ask for a tiling interpretation for the weighted FM-paths, describing Ra.n as a function of the seed data 

XM- 

Consider first Mmax, with ma = a — 1, where i?i.„ is a partition function of weighted Dyck paths on 
the strip < y < 2r +1 from (0, 0) to (2n, 0). Comparing this with the situation described in Lemma [9731 
we can make a direct connection by forbiding the steps (2, 0), or by using tiles a, b, c only. Thus, Ri^n/Ri,o 
is the generating function for tilings of the domain HAn by means of tiles a, b, c only, and with weights 1 
per tiles a or c, and yi per tile b with center at height i — 1/2, for i — 1,2, 2r + 1. 

More generally, let Fm — 72r+2(*ii •■•j is): a skeleton tree as in Definition 15.71 As before we choose 
the weights Wq — wi ^ 0, Zj = for j > 2{2r -I- 1 — s) — 1, Z2(2r-s) = 2/2r+i(M), and Z2j-i = for all 
indices j G [2, 2r — s — 1] \ {ii, 12, is}, corresponding to missing horizontal edges at height j in Fm- The 
remaining z^'s are identified with yj(M), in increasing order of indices. Again, the vanishing conditions 
impose some truncation to a strip h E [—1/2, 2r — s + 3/2], and forbid the use of tiles d except at heights 
ii,«2, 

Next, consider the case of the Motzkin path M(fc) := /ii/^2---/^fc(Mo), with mi — TO2 = • • • = = 1, 
and mj = for j > fc. The graph FM(fe) is Gr with one additional down-pointing edge fc + 2 — > fc, assuming 
that r > fc + 1. 

We define modified Aztec paths (A]y[(-fe-)-paths) as the corresponding F]y[(-fc)-paths in the quadrant 
(without the usual restriction due to r being finite). It has an extra step (0, —2) only from height h = k + 2 
to height h = k. The weighting is as in Section and the extra step receives a weight Zk+2,k- Paths 
from (0, 0) to {2n, 0) are in bijection with tilings of the domain depicted in Fig|32l by means of the domino 
tiles of Fig|28l plus rigid pairs of 1 x 1 defects tiles (see the medallion in Fig|32|), with centers at positions 
{t — 1/2, k) and {t + 1/2, k + 2) with k + 2<t<2n — k,t integer. The "maximal" path pq now involves 
a succession of n + 1 up steps (1, 1), followed by n — fc down steps (1, —1), then one down step (0, —2) 
and finally fc down steps (1, —1). The local moves (i) and (ii) of Fig |30l must now be supplemented by new 
moves expressing the reduction of the special down steps (0, —2). 

The partition function for weighted AM(fc)-paths is identified with i?i^„ expressed as a function of XM(fc) 
(up to a multiplicative factor of i?i,o), choosing the weights wq — wi — Z2k+i — Z2k+3 — 0, Zj = for 
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2n-2 



Figure 32. The bijection between modified Aztec j4M(fc)-paths from (0,0) to (2n, 0) and the 
domino tilings of a domain made of the half Aztec diamond (dashed lines), enhanced as indicated, 
tiled by means of dominos and rigid pairs of 1 x 1 squares with centers at height k + 2 and k 
(depicted in the medallion). 
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Figure 33. The local rules for constructing the domain T)^ out of the maximal j4M-path. 

j > 2r + 4, Z2fe+2 = J/2fc+i, Zj = Vj for j = l,2,...,2fc, and Zk+2,k = Vk+2,k = V2kV2k+2/V2k+i, where 
the Uj = yj(M(fc)). Equivalently, Ri^n/Ri,o is the partition function for weighted tilings of the domain 
of Figure [32l delimited by the tiles corresponding to the maximal path , and by the floor at height 
h = -1/2. 

For an arbitrary Motzkin path M, the Am Aztec paths are the Fivi-paths without restriction due to r. 

(n) 

Definition 9.5. The domain is defined as follows. Starting from Fm, remove any down-pointing 

edges i ^ j */ Fm ho-s an edge i' — > j' with [j,i] C (strict inclusion). This leaves us with a set of 

"maximal" down-pointing edges. Define the "maximal" path from (0,0) to (2n, 0), which has no horizontal 
step, and goes as far up and to the right as possible, namely starts with a maximal number of steps (1, 1), 
descends via n maximal descending steps. The domain is constructed by associating boundary pieces to 
each step, as shown in Fiaure lSM Finally, the domain is completed by a horizontal line at height h — —1/2. 

Lemma 9.6. The AM.-po,ths from (0,0) to (2n, 0) are in bijection with tilings of the domain by 
means of the usual 2x1 and 1x2 tiles, plus rigid pairs of square 1x1 tiles centered at points of the form 
(t — and (t + ^—^^,i) for all the pairs i,j of vertices connected by down-pointing edges on Fm- 

We refer to Figure |34] for an example. We can now say that Ri,n+mi/ Ri.mi is the partition function of 
tilings of 2)^^ as described above. 

Finally, recall that Ra,n+mi / {Ri,mi )" is the generating function for families of a strongly non-intersecting 
Fjvi-paths, starting at (0, 0), (2a — 2, 0) and ending at (2n, 0), (2n + 2a — 2, 0). These paths can be 
represented as tilings of the domain with a — 1 additional indentations on the bottom left and bottom 
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Figure 34. A typical domain , with its maximal Ajvi-path (in red), the target graph Fm 
and the Motzkin path M. We have indicated the two maximal edges, along which the maximal 
path takes its maximal descents. 




Figure 35. The indented domain 1)^ , serving for the representation of Ra,n in terms of tilings, 
and the target graph Fm (here a — 3). Note that the tiling involves defect pairs corresponding 
to descents of Fm, and that horizontal edges (tiles d) are forbidden at heights where Fm has no 
horizontal edges. Outside of the shaded blue domain, all the tiles are fixed by the indentations 
and the absence of d tile in the bottom row. 

right boundary (see Fig|35p. by means of dominos and rigid defect pairs corresponding to down-pointing 
edges of Fm- The strong non- intersection of the paths imposes extra constraints on the tilings, by forbiding 
some local configurations. 

10. Conclusion 

In this paper, we found a simple structural explanation for the cluster mutations, which allow for 
sweeping the set of possible initial data for the Q-system, in terms of simple local rearrangements of 
the continued fractions that generate the i?i^„'s. This local move allowing for generating mutations is 
reminiscent of the local "Yang-Baxter" -like relation used in [8] in the context of total positivity of the 
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Grassmannian, as expressed through local positive transfer matrices for networks. We do not yet fully 
undertand this relation, although a partial explanation is found in [5 . 

In view of Example 16.81 we can think of the constructions of this paper as generalizations of the Stieltjes 
theorem [53], that now allow to re-express the series F{X) = X]fc>o(^-'^)'''^'=/'^'^^^ different ways as 
(mutated) possibly multiply branching continued fractions, whose coefhcients are particular combinations 
of Hankel determinants involving the sequence a^, each such rewriting corresponding to a Motzkin path. 
To make the contact with our results, we simply have to take t = — 1/A, Ofe — Ri^k and to identify 
Ra,n = AJ^~"+^ with the Hankel determinants of eq. (j6.27|) . For each Motzkin path M, we find a new 
continued fraction expression for F{X) involving only the Hankel determinants corresponding to the cluster 
variable xm, via the weights j/i(M). 

Equations called Q-systems exist for all simple Lie algebras. We have checked that in these cases, 
also satisfies linear recursion relations with constant coefficients. We expect the constructions of the present 
paper to generalize to all these cases. In particular, we expect cluster positivity to be a consequence of the 
LGV formula applied to counting possibly interacting families of non-intersecting paths. 

The Q-system is a specialization of the T-sytem, a discrete integrable system with one additional pa- 
rameter. It was shown in [3] that T-systems can be considered as cluster algebras, in general, of infinite 
rank. The statistical models in this paper are particularly well-suited to the solution of the T-system [B] . 

In the case of the cluster variables corresponding to seeds of the cluster algebra outside the subgraph 
9r, we do not have an interpretation in terms of a statistical model. Outside this subgraph, the exchange 
matrix B has entries which grow arbitrarily. This seems to suggest that the evolution in directions leading 
outside Sr is not integrable. 

Appendix A. The case of A3 in the heap formulation. 

In this appendix, we detail the heap interpretation of the solution Ri_„ of the ^3 Q-system, as expressed 
in terms of the various cluster variables corresponding to the 9 seeds of the fundamental domain. We 
apply the same systematic method as in Example 14.81 For simplicity, we have just depicted in FiglHSl 
the heap graphs, together with the mutations relating them, for each cluster variable of the fundamental 
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domain. The corresponding rearrangements of the generating function {t) and new weights are given 
below, with the labehng of vertices indicated in Figl36l These are straightforwardly generated by carefully 
following the successive applications of Lemmas 14.51 and 14.61 

Initial cluster variable: G3. We start from the graph G3 and the associated heap generating function 
F^^\t) = 1 + with 



with the y's as in p.l7p . namely 

_ Ri,i _ R2,i _ fii, 0^2,1 _ ^1,0^3,1 _ R2,qR3S _ R2,o _ Ra.Q 
Ui — 15 I 2/2 — ^5 ,2/3 — ^5 : J/4 — ^5 ,2/5 — -f^ ,2/6 — -f^ J5 ,2/7 — -f^ ■ 

J^l.O J^ifiUi,! 0^1,1 -rt2, 0^2,1 -K3,o-K2,l ^3,0-^3,1 K3_i 

Note that the weights are related via yiyzyhUr ~ 1- 

Mutation iii{G^). We apply Lemma [4.61 to vertices 1, 2 and the structure attached to 3. This gives 
F[^\t) = 1 + with 



1-t- 



^ 1 — t27 

with 

2/22/3 2/22/4 2/12/3 2/12/4 

^1 = 2/1 + 2/2, 22 = ■ ,28 = \ , 2^3 = \ ,2^4 = ^ ,2^5 = 2/5, 2^6 = 2/6, ^^T = 2/7 

2/1 + 2/2 2/1 + 2/2 2/1 + 2/2 2/1 + 2/2 

namely, upon using the Q-system relation 0^1,2 = ^11+ ^2,1: 

-R1.2 -^2,1 ^1.1-^2,1 ^1,1-^3,1 



Z-\ = — , Zo — ■ , Z3 = ■ ,Za — 

JrC2,0-rii.iKi.2 ^2,0-^1,2 ^2,0^2,l-fll,2 

^2,0^3,1 ^2,0 ^3,0 ^3,1 

^5 = , 26 t; , Z7 = — , Zs — 



-R3, 0^2,1 -R3, 0^3,1 ^3,1 R2flRl,2 

Note that there is one more z than ?/'s, but that we now have two relations 2:1 03^5^7 = 1 and z^zg, = ziz/^. 
Mutation [i2{G-i). We apply Lemma ITBl to vertices 3, 4 and the descendent structure attached to 4. This 
gives Ff'it) = 1 + t?/i$(^2(<^3))(t)^ with 

(M2(G3))/.N _ \ 



1 — tx\ — t- 



with 



2/42/5 2/42/6 2/32/5 2/32/6 

a;i = yi,X2 = y2,X3 = 2/3 + 2/4,2:4 = ■ ,a;8 = ■ ,2:5 = ^ ,xe = ■ ,^7 = yj. 

2/3 + 2/4 2/3 + 2/4 2/3 + 2/4 2/3 + 2/4 

namely, upon using the Q-system relation i?2,o^2,2 — R2 1 + Ri,iR3A- 

Ri,i R2.1 Ri,qR2,2 Ri^Rl.i 



-Ki,o -Ki,o-rti,i -n-i,irt2,l -rt3 0^2,1/12,2 

R2,lR'i.l R2A Rsfi Rl,l 
^5 = 7^ ;7~,2;6 = -fi j:; ,XJ = — ,Xs — 



/C3, 0/1-2, 2 /<-3, 0/13, 1/1-2, 2 /t3,l /i-3,o/t,2,2 

As for the z's above, there is one more x than y's, but they obey two relations: xix^x^xj ~ 1 and 



X5XS = a::4a:6- 
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Mutation fi^{G3). This is part of the sequence of graphs of Example 15.61 It is obtained by applying the 
Lemma IT^ to the vertices 5, 6, 7 of G3. This gives F^^\t) = 1 + ij/i^^^^lGa)) ^^^^^ 



i-tti~t- 



with 



2/62/7 2/52/7 
tl = 2/1:^2 = 2/2,^3 = 2/3,^4 = 2/4,^5 = 2/5 + 2/6-^6 = \ ,h 



2/5 +2/6 2/5 + 2/6 

namely, upon using the Q-system relation R^.qRz,2 = ^31+ R2,i'- 

, _ Ri.i _ R2,i _ Ri,oR2,i _ RifiRs.i , _ -R2, 0-^3,2 _ R2,i _ Ra.i 

^ ~ ~Q ' ^ ^ 1? P ' ^ "r R ' ~ "p R ' ^ ^ "r P ' ^ ~ "p R ' ~ "r ■ 
^1,0 -ttl.O-fl-ia -fl2,0^ia -n-2, 0^2,1 ^2,1-^34 -tt3^iK3_2 ^3,2 

Note that the t satisfy the relation tit^t^ti — 1. 

Mutation /ii^2(G3). We apply Lemma to the vertices 1, 2, and the substructure attached to 3 in the 
graph fi2{G3) (see Figl36l)- This gives $(A'iA'2(G3))(i) = $(A'2(G3))(i)^ ^[^^^^ 

(A.l) $(M1M2(G3))(^) _ 1 



1 + "1 



"4 + 1-^^ 



with 



a;2a;3 2:10:3 

Ui = Xi + X2, U2 = , M3 = , U4 — X4, Ms = 2:5, Uq = Xq, M7 = Xj, Us — Xs- 

Xi + X2 Xi+ X2 



namely, upon using the Q-system relation 0^1,2 — Rf 1 + ^2,1 



i?1.2 i?2,2 Rl,lR2.2 -Rl, 1-^3,1 



Ui — —.U2 — , U3 — — ,W4 

P' PR' PR' RPR 

-Tlia J^isni,2 -rtl 2-K2,l -fl3,o-ri2,lK2,2 

-R2, 1^3,1 ^2,1 -^3,0 -Rl.l 

"5 = 7^ n ,""6 = 7:; ,U7 = 77— ,"8 



Rz,oR2a ' -R3, 0-^3,1-^2, 2 ' -^3,1 ' -R3, 0-^2,2 

Note that the u's satisfy the two relations uiu^u^u^ — 1 and itsug = u^ue. It is instructive to see how to 
arrive at the same result from the sequence of mutations /i2/ii(G'3) = /^i/i2(G3). We now start from the 
graph fii^Gs), and apply the Lemma 14.61 with (i) the vertex 2 and its substructure, (ii) the vertex 8 and 
its substructure common with that of vertex 2 (iii) its attached substructure not common with that of the 
vertex 2, namely with 

Z2 



1 - tZ3 - tj 



Z8 



C = 



1 - tz^ 



The result is 

(A.2) $(M2Pi(G3))(^) 



\-t- 
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with 



Z2 + 28 
2:22:5 



y2 = Ml, -"2 



22 + Z8 = J/2 



XT-, 



y4 



Z2Z6 



yi + y2 

Z9,Z(, _ 
Z2 + Z8 

ysye 



X2X3 

Xi + X2 
2/42/6 



U2, 



2/3 



2/4 



xs = Us, 



zj = xr = M7, 



Z2 + ^8 Z2 + zs ys + 2/4 

where we have used the expressions of the z's and the x's in terms of the y's. Applying again Lemma 14.61 

Z3, 6 = Z4, c = Z5 + zq/{1 — tzi), we may rewrite: 



to the first factor under u'2, with a 



(A.3) 



with 



23 + 



Zi 



1 — — 



1 - t- 



23 + 24 
2325 



2/1- 



2/3 + 2/4 

2/1+2/2 

2225 



W5 = U5,W7 



Z4Z5 



2325 



23 + 24 22 + 28 



W4 = W4, 



27 = M7, 



23 + 24 22 + 28 

where we have used the relation 23/24 = 22/23 between the 2's. Substituting this into (|A.2p . we finally 
recover eq. ljA.ip . This case is made particularly cumbersome by the fact that the vertex 2 contains a 
substructure (which has not yet been the case so far), encoded in the fraction (jA.3|) . From this example, 
we learn that the mutation ^2 consists actually of two simultaneous applications of the Lemma 14.61 one 
"usual" and one within the substructure attached to the initial vertex 2. 

Mutation 111(13(03). We apply Lemma l46l to the vertices 1, 2, and the substructure attached to 2 in the 
graph ^3(03) (see FiglMl)- This gives ^(A'lMslGa)) = $(^'3(63)) (^)^ ^[^]^ 

(A.4) $(mim(G3))(^) 1 



l~t- 



"2 + - 



with 



Wi 



h + t2,W2 



t2t3 



^2^4 



,W3 



tlt3 



tl+t2 ti+t2 ti+ t2 

namely, upon using the Q-system relation R1.0R1.2 = Ri i + R 



tl+t2 



Wi 



W5 



-,W2 = 



Ri^ 

R2fiRs,,2 

R2AR2,, 



R2.1 



R: 



R2,oRi,iRi,2 



W3 



3,1 



R 



2,1 



-,W7 



Ri,2R2,a 
R3.1 



,W4 



Rl 1-^3,1 



R2,oR2,lRl,2 



^3,1-^3,2 -^3,2 R2,qRi,2 

Note the relations: wiW3W5W'j = 1 and W2Wa = 'W3Ws- For completeness, let us compute the function 
$(ai3/^i(G3)) ^ $(/'i(G3))^ a,nd check that indeed ^1/^3(03) = /X3Ati(G3) yield the same result. We now apply 
the Lemma l46l with the vertices 5, 6 and 7 of the graph /ii(G3), with respectively a = 25, b — zq and 
c = 27, while d — 0. This gives 

(A.5) $(^3/^1(63)) ^ I . 



1-t- 



n' 4- _ 
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with 

, , 1/21/3 ^2^3 

Wi = Zi ^ yi + y2 = ti + t2 = Wi, — = - — — - = W2 

yi +1/2 ti+ t2 

I V2Vi t2U I 2/12/3 hh 

'Ws = Zs^ ■ = — — - = Ws, ^3 = Z3 = ; = — — - = W3, 

2/1 +2/2 h+ h 2/1 +2/2 ^1 + h 

I 2/12/4 ^1^4 / 
W4 = Z4 = ■ = — - = Wi, W5 = ^5 + Z6 = 2/5 + 2/6 = ^5 = W5, 

2/1 +2/2 ti + h 
I Z(iZ-j y^y-j , Z5Z7 yzy^ 
Wq = ■ = ■ = ^6 = We, w-j = ■ = ■ = = u't, 

^5 + Z6 2/5 + 2/6 ^^5 + 26 2/5 + 2/6 

where the various identifications are made by use of the expressions of the z's and i's in terms of the y's. 
We deduce that the expressions (|A.4I) and (|A.5I) are identical. 

Mutation /Z2//3(G3). This is part of the sequence of graphs of Example 15.61 It is obtained by applying 
the Lemma l4!6l to the vertices 3,4 and the structure attached to 5 of the graph /i3(G'3) (see FiglM]). This 
gives $('^2P3(G3))(^-) = $(m3(G3))(^)^ ^-^^-^ 

(A.6) $(M2A.3(G3))(^)^ \ ^ 

with 

^4^5 ^3^5 

Sl =il,S2 =^2,53 =^3+^4,54 = " — r'*5 = T — ^'■^6 =^6,57 = . 

H +14 13 + 14 

namely, upon using the Q-system relation R2.qR2,2 = ^21+ ^i,i^3.i- 

_ -Rl,l _ R2,l _ RlflR2,2 _ -Rl,1^3,2 _ R2,lRz,2 _ R2,l _ Rj,! 

~ ~P ' ~ ~R R ' ~ ~R R ' '^^ ~ "r R ' ~ "r R ' '^^ ~ "r R ' ~ 1? ' 
-n-1,0 ^1,0^1-1,1 ^1,1-^24 ^2,1^12,2 ^3,1^2,2 ^3,1JL3^2 -^3^2 

Note the relation S1S3S5S7 = 1. For completeness, let us compute the function $('^3''2(G3)) _ <j)(M2(G3))^ 
and check that indeed H2f^3{G3) — ^^3^2(^3) yield the same result. We now apply the Lemma on the 
graph /i2(G3), with (i) the vertex 4 and its substructure (ii) the vertex 8 and its substructure common 
to vertex 4, and the substructure of vertex 8 not common to vertex 4. This amounts to taking a = 
X4/{1 — tx5 ~ txQ/{l — tx-j)), b — xs/{l — tx5 — txe/{l — tx-j)) and c = x-j, while d = 0. This yields 

1 



<J)(A'3A'2(G3)) 



1 - is'i - t- 



i~t- 



with 



s'l = xi ^ yi =ti ^ Sl, s'2 ^ X2 = y2 = h = ^2, S3 = 2:3 = 2/3 + 2/4 = <3 + ^4 = S3, 

, , 2/5 + 2/6 ^4*5 , 2:7X8 yayr 
S4 = a;4 + = 2/4 : = t — — - = S4, Sg = ■ = ■ = ^6 = sq, 

2/3 + 2/4 + ti X4 + Xs 2/5 + 2/6 

, X7a;4 2/52/7 

tj = S7, 



Xi + XS 2/5 + 2/6 

where identifications are made by expressing the x's, s's and t's in terms of the y's. As in the case of 
^2^1(^*3), we see that we must still apply the Lemma 14.61 to the first factor under s^ in (jA.Tp . which 
corresponds to the substructure attached to the initial vertex 4. This allows to rewrite 



X5 

with 



Xq 



1 — tXv 



2/5 + 2/6 ^3*5 „ XeX7 yeyv 

•.5+Xq = 2/3 ■ = 7 — - = S5, Sg = ■ = ■ = tg = S6 

2/3 + 2/4 ^3 + <4 2:5 +Xe 2/5 + 2/6 

a;5a;7 2/52/7 , 

= ^7 = S7- 



2^5 +2:6 2/5 + 2/6 
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Substituting these into (|A.7[) allows to identify it with (jA.6[) . 

Mutation ^6/^2/^3(^3) = Hi,. This is part of the sequence of graphs of Example 15.61 It is obtained by ap- 
plying the Lemnia|4i6]to the vertices 5, 6, 7 of the graph /X2/i3(G3) (seeFiglSSl). This gives ^(^'^''^^(Ga))^^-) ^ 
$(/x2/^3(G3))(i)^ with 



l-/ri-/— 



with 



ri 



si,r2 = S2,r3 = S3,r4 = S4,r5 = 55 + 56,^*6 



S6 



S5S7 



S6 



namely, upon using the Q-system relation i?3.i^3,3 — R3 



i?2 



ri 



,r2 



R2A 



RifiRi.i 



Ri,qR2.: 
R1SR2, 



RiaR. 



3,2 



-^2,1^2,2 



R2,lR: 



■3,3 



-R2, 2-^3,2 



,r6 



R2 



^3,2-^3,3 



,r7 



R: 



3,2 



R 



3,3 



Note the relation rir^r^r^ — 1. 

Mutation /i4/ii/t2(G3). It is obtained by applying the Lemma [4.61 to the vertices 1,2 and the structure 
attached to 2 of the graph ^1(12(03) (see FiglMl). Writing $(^'iA'2(G3))(^) = 1 + tui^'^^i^^C^^)) (i)^ .^^e have 

XJ,(M4A'1A'2(G3))('^^ _ XI|(M1M2(G3))^^-)^ .^^j^J^ 



(A.8) 



^(/J4/J1A'2(G3)) ^^-J 



1 - <- 



1-t 



with 



U2U3 

Wl = Ul + U2, W2 = ; , V9 



U1+U2' 



U2U4 
U1+U2' 



V4, 



U1U4 

Ul + U2 



, Wio = ^ ,W3 

Ul + U2 



MlM3 
Ul + M2 ' 



Ml + U2 



namely, upon using the Q-system relation i?i,ii?i,3 = i?i 2 + ^2,2: 



,V2 



Rla 



ve = 



Ri^ 

Ri,2 ^2,1^1,2^1,3 

-^2,1 



,V3 



^1,2^2,2 



-R3, 0-^3,1^2.: 



^3.0 

V7 = — — ,Vs = 
-K3.1 



R2,lRl,3 
,2 



P2 r2 
^1,2-^34 



-R3. 0-^2,1^2,2^1,3 

4 



^3,0^2,2^1, 



-,Vg 



RIa 



R3.oR2,lRl.: 



,V5 



,Vl0 = 



R2,lR3,l 
^3,0^2,2 
1 

-R3.0-Rl,3 



Note that as we have 3 more w's than y's, we expect four relations between them. These read: W1W3W5W7 — 1, 
V2VA = W3W9, = W4W10 and i;2W8 = ^swio- 

We conclude that for r = 3, i?i,n is a positive Laurent polynomial of all the mutations of the initial data 
xo in the fundamental domain 3^3. 



Appendix B. The case of ^3 in the path formulation. 

In this appendix, we detail the path interpretation of the solution Ri^n of the A3 Q-system, as expressed 
in terms of the various initial data corresponding to the 9 cluster variables of the fundamental domain J's, 
labelled by the 9 Motzkin paths of FigH 

As a direct illustration of the constructions of Section |6l we list below for each of the nine cluster 
variables the corresponding transfer matrix Tp, and the generating function ((/ — rr)^^)g qi in terms of 
dummy weights yi and yiJ^ that are actually shorthand notations for the weights 2/i(M) and their redundant 
counterparts, to be extracted from Theorem 16.71 and eq. (|5.5p . The corresponding graphs are represented 
in Figi37l together with the mutations relating the associated cluster variables. 
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For Mq = {(0, 1), (0, 2), (0, 3)}: The corresponding target graph is Fmo = ^3. The transfer matrix and 
the associated generating function are: 
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/ 



Mq) 



1 - t 



i-f- 
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For Ml = /xiMo = {(1, 1), (0, 2), (0, 3)}, with yg^i = 2/2^4/2/3: 
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tyi 
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0/ 



i-t- 



V2 + 



V3,l 
l-t!/5-*T 



-ty2 - 



ye 
-tyj 



tV7 



For M2 = /X2M0 = {(0, 1), (1, 2), (0, 3)}, with ^5,3 = 2/4j/6/j/5: 



'Pm, 



'Pm, 



/O tyi 
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\0 
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MsMo = 


{(0,1), (0,2), (1,3)}: 




/O tyr 
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35: 



For M2,i = M2M1M0 = {(1, 1), (1, 2), (0, 3)}: With ^4,2 = 2/42/6/2/5, 



- M2,l 
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For M3,i = /X3/U1M0 = {(1, 1), (0, 2), (1, 3)}: With 2/3,1 = 2/2^4/2/3, 
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For M3,2 = A'3A'2Mo - {(0, 1), (1, 2), (1, 3)}: 
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For M6,3,2 = Ai6Ai3Ai2Mo = {(0, 1), (1, 2), (2, 3)}: 
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Note finally that the graphs in Figures l36l and [37l are duals of each other. This duality is best seen by 
expressing a bijection between the paths (represented with their time extension) and the corresponding 
heaps, by associating a disc of the heap to each descent of the path. The heap graph is simply the graph 
whose vertices stand for the discs and whose edges indicate the pairs of overlapping discs. 
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